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The dynamical breaking of chiral symmetry in QCD is caused by nonperturbative interactions on 
a distance scale p ~ 0.3 fm, much smaller than the typical hadronic size R ~ 1 fm. These short- 
distance interactions influence the intrinsic transverse momentum distributions of partons and their 
correlations at a low normalization point. We study this phenomenon in an effective description of 
the low-energy dynamics in terms of chiral constituent quark degrees of freedom, which refers to the 
large- N c limit of QCD. The nucleon is obtained as a system of constituent quarks and antiquarks 
moving in a self-consistent classical chiral field (relativistic mean-field approximation, or chiral 
quark-soliton model) . The calculated transverse momentum distributions of constituent quarks and 
antiquarks are matched with QCD quarks, antiquarks and gluons at the chiral symmetry-breaking 
scale p~ 2 . We find that the transverse momentum distribution of valence quarks is localized at p% ~ 
R- 2 and roughly of Gaussian shape. The distribution of unpolarized sea quarks exhibits a would- 
be power-like tail ~ 1/Pt extending up to the chiral symmetry-breaking scale. Similar behavior is 
observed in the flavor-nonsinglet polarized sea. The high-momentum tails are the result of short- 
range correlations between sea quarks in the nucleon's light-cone wave function, which are analogous 
to short-range NN correlations in nuclei. We show that the nucleon's light-cone wave function 
contains correlated pairs of transverse size p <C R with scalar-isoscalar (E) and pseudoscalar- 
isovector (II) quantum numbers, whose internal wave functions have a distinctive spin structure 
and become identical at p% ~ p~ 2 (restoration of chiral symmetry). These features are model- 
independent and represent an effect of dynamical chiral symmetry breaking on the nucleon's partonic 
structure. Our results have numerous implications for the transverse momentum distributions of 
particles produced in hard scattering processes. Under certain conditions the nonperturbative parton 
correlations predicted here could be observed in particle correlations between the current and target 
fragmentation regions of deep-inelastic scattering. 
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I. INTRODUCTION 

The parton model provides the basic script for express- 
ing hadron structure as seen by short-distance probes 
such as local current operators or high-momentum trans- 
fer processes. Its fundamental assumption is that the 
fast-moving hadron can be regarded as a collection of 
pointlike constituents that behave like free particles on 
the timescale of their interaction with the external probe. 
The basic quantities are the number densities of partons 
as functions of their longitudinal momentum fraction x. 
Matrix elements of local current operators measure in- 
tegrals of the parton densities ("sum rules"), while pro- 
cesses such as deep-inelastic eN scattering (or DIS) or 
production of high-mass systems in NN scattering probe 
them differentially in x. 

The parton model can be thought of as the limiting 
case of a large class of dynamical models of hadron struc- 
ture, in which the typical transverse momenta do not 
grow (or at most logarithmically) as the longitudinal mo- 
mentum of the hadron is increased. In this approach the 
parton densities appear as integrals over the transverse 
momenta of the constituents [l|, 

fi(x\f) = [ d 2 PT fi(x, PT ), (1.1) 

where /i 2 signifies a cutoff (e.g. in the parton virtuality, 
or in the invariant mass of configurations in the light- 
front wave functions), which restricts the integral over 
Pt and defines the resolution scale at which the picture 
of pointlike partons is supposed to apply. The integrand 
in Eq. is referred to as the "intrinsic" transverse 

momentum distribution of the partons. It is tempting 
to interpret this function as a density of particles and 
use it to study the internal motion of the constituents in 
the system. While this is principally possible, studies of 
field-theoretical models show that the transverse motion 
can generally not be separated from the interactions in 
the system. In gauge theories such as QED and QCD the 
parton picture requires the light-cone gauge, where the 
transverse components of the gauge potential, At, repre- 
sent dynamical degrees of freedom. Gauge invariance im- 
plies that the transverse derivatives of the fermion fields 
appear in the combination Vt = 8t — iAt, linking the 
kinetic transverse momentum of the fermions to dynam- 
ical gauge fields in the hadron 0]. The interpretation of 
the intrinsic transverse momentum distributions is thus 



generally much more subtle than that of the parton num- 
ber densities. 

A partonic picture of hadron structure is expected to 
emerge from QCD as the result of long-range nonpertur- 
bative interactions, and is commonly used to describe the 
boundary conditions for pcrturbative QCD calculations 
of hard processes. In this approach the parton densities 
can be expressed as matrix elements of certain quark and 
gluon light-ray operators of twist 2, i.e., correlation func- 
tions of the fields at light -like separation, normalized at 
the scale fi 2 (normalization point) Q. The dependence 
on the scale can be calculated pcrturbativcly for suffi- 
ciently large values; it is governed by the renormaliza- 
tion group equations for the composite operators, which 
coincide with the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (or DGLAP) evolution equations describing parton 
decay in the leading logarithmic approximation. In the 
region where perturbative evolution is applicable the in- 
trinsic pt distribution of partons is well-defined and of 
the form 

h{x,pr) ~ (1.2) 
P T 

which reflects the ultraviolet (or UV) divergences of QCD 
and implies a logarithmic dependence of the parton den- 
sity Eq. (|l.ip on the scale /j, 2 . Conversely, the coefficient 
C/j (x) can be recovered as the logarithmic derivative of 
the parton density with respect to the scale, 

C h {x) = tt-V^j fMf)- (1-3) 

"Unintegrated" parton densities defined according to 
Eqs. (|1.2p and (|1.3|) have been employed in phenomeno- 
logical studies of DIS at small x (< 10 -2 ) j4[. A more 
careful treatment accounts also for the loss of partons in 
the course of DGLAP evolution as described by the Su- 
dakov survival factor, which becomes more important at 
larger values of a; [H, H . 

The extension of the concept of parton transverse mo- 
mentum in QCD into the nonperturbativc domain is in- 
herently not unique, and different prescriptions have been 
proposed in the context of studies of various classes of 
high-momentum transfer processes. The calculation of 
power corrections to unpolarizcd DIS structure functions 
in the collinear expansion leads to twist-4 light-ray op- 
erators of the form ip^T,i ■ ■ ■ ^Tjip, which can be in- 
terpreted as measuring the "average transverse momen- 
tum" of the quarks (keeping in mind that it cannot be 
separated from the transverse gauge fields) 0, 01 ■ Sim- 
ilar operators of twist 3 appear in the study of single- 
spin asymmetries in high-pr particle production Q. An 
alternative prescription are the transverse-momentum- 
dependent distributions (or TMDs) introduced in the 
context of a QCD description of semi-inclusive DIS at 
low transverse momenta [9l4llj . They are defined as cor- 
relation functions of fields off the light-cone, at finite 
transverse separations, and involve gauge links (or phase 
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factors) describing the effect of QCD initial/final-state 
interactions of the quark participating in the hard pro- 
cess. The renormalization properties of these operators 
and the proper choice of gauge links are presently the 
subject of intense theoretical study (l2l - fl4( . The formal 
correspondence between the "collinear" and "TMD" def- 
initions of parton transverse momentum in the perturba- 
tive region of high pt was studied in Ref. fl5l | . A natural 
concept of "intrinsic" transverse momentum also appears 
in approaches which describe hadron structure in terms 
of light-cone wave functions at a low scale 0, [I?} . 

The effective dynamics of strong interactions at the 
hadronic scale is in large measure governed by the spon- 
taneous breaking of chiral symmetry in the QCD vac- 
uum. It leads to the appearance of a chiral conden- 
sate (tpip) (order parameter), whose phase fluctua- 
tions give rise to almost massless excitations (Goldstonc 
bosons), the pions. Their interactions are summarized 
by the universal chiral Lagrangian and determine the be- 
havior of strong interactions over distances of the order 
~ 1/tojt, which can be studied using methods of effective 
field theory. However, the influence of chiral symmetry 
breaking extends down to much shorter distances, defined 
by the range of the nonperturbative QCD interactions 
that lead to the appearance of the chiral condensate. Lat- 
tice QCD calculations and numerous phenomenological 
observations suggest that the size of the nonperturbative 
field configurations causing the spontaneous breaking of 
chiral symmetry in QCD is much smaller than the typical 
hadronic radius R ~ 1 fm. An objective gauge-invariant 
measure of this scale is the average quark virtuality in 
the chiral condensate [HI 

2 (^) ■ { ■ > 

Lattice simulations give mg/2 > 0.5 GcV 2 at a normal- 
ization point of (j. ~ 1 GeV [l9|, [2(| ; even larger values 
were obtained in Ref. [2l|. This is supported by the re- 
sults of direct studies of chirality-ffipping topologi cally 
charged vacuum fluctuations in lattice simulations [22j; 
see Refs. 0, HH for a review. The same conclusion is 
obtained from the instanton model of the QCD vacuum 
IH [Hi (see Refs. [H, H3| for a review), where the typi- 
cal size of the instantons is p ~ 0.3 fm and the average 
quark virtuality in the chiral condensate is found to be 
ml/2 = 2p~ 2 « 0.7 GeV 2 [Hj]. Abstracting from these 
findings, one may state that the spontaneous breaking of 
chiral symmetry in QCD is characterized by a dynamical 
scale much shorter than the typical hadronic radius, 

p ~ 0.3 fm < R ~ lfm. (1.5) 

The existence of this nonperturbative short-distance 
scale has far-reaching consequences for the structure of 
hadrons and their low-energy interactions [H H3] ■ 

Here we want to ask what dynamical chiral symmetry 
breaking and the existence of the short-distance scale 
p imply for the transverse momentum distribution of 



partons in the nucleon at a low normalization point. 
This question is clearly of great importance for both the 
general theoretical understanding of partonic structure 
and the phenomenology of hard processes with identified 
particles, such as semi-inclusive ep scattering, jets and 
Drcll-Yan pair production in pp scattering, and multipar- 
ton interactions in pp collisions. In view of the ambigui- 
ties in the very definition of intrinsic transverse momen- 
tum in QCD we shall not attempt to approach this prob- 
lem in a model-independent manner, as the evaluation 
of certain a priori defined QCD operators in the nucleon 
state. Instead, we shall study the transverse momentum 
distributions of partons in a model of the effective low- 
energy dynamics resulting from the spontaneous breaking 
of chiral symmetry, which implements the two dynamical 
scales of Eq. (|1.5I) . The dynamical model will suggest a 
natural definition of the intrinsic transverse momentum 
distribution in terms of effective degrees of freedom, in- 
cluding the pertinent resolution scale. The matching of 
the model pr distributions with QCD will then be con- 
sidered on the basis of their specific form, and with the 
help of empirical information on the p^-integratcd par- 
ton densities, at a normalization point determined by the 
chiral symmetry-breaking scale, p 2 ~ p 2 - In the present 
state of development such an approach is fully justified 
and provides a useful complement to more abstract stud- 
ies of transverse momentum distributions based on spe- 
cific QCD operator definitions. 

Numerous observations point to the importance of con- 
stituent quarks and pions as effective degrees of free- 
dom below the chiral symmetry-breaking scale. Theo- 
retical arguments suggest that in the large-_/V c limit of 
QCD the effective dynamics resulting from the sponta- 
neous breaking of chiral symmetry can be approximated 
by a field-theoretical model based on chiral constituent 
quarks [H [2(| . It expresses the fact that the modes of 
the QCD quark fields with virtualities below the chiral- 
symmetry breaking scale p~ 2 acquire a dynamical mass. 
Because of chiral invariance, this is necessarily accom- 
panied by a coupling to the Goldstone pion field, which 
in the large-iV c limit is itself a composite of constituent 
quarks and antiquarks. This effective dynamics is rel- 
evant up to the chiral symmetry-breaking scale, which 
appears as the UV cutoff of the model. A crucial point is 
that the chiral symmetry-breaking scale is assumed to be 
parametrically large compared to the dynamical quark 
mass, such that the massive constituent quarks can be 
regarded as pointlike over a wide range of virtualities. 
This two-scale picture gives a precise meaning to the no- 
tion of constituent quarks as effective degrees of freedom 
and provides an ordering principle for the calculation of 
hadron structure. 

In the effective chiral model the nucleon is obtained 
as an extended solution in which massive quarks and 
antiquarks move in the background of a self-consistent 
pion field (relativistic mean-field approximation, or chi- 
ral quark-soliton model) [30| . Matrix elements of oper- 
ators between nucleon states can be computed in a sys- 
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tematic 1/N C expansion. This picture results in an es- 
sentially parameter-free description of the static nucleon 
observables and form factors [31(. Because the descrip- 
tion is fully field-theoretical, the model has a partonic 
limit and can be used to calculate the nucleon's parton 
densities at a low normalization point [32l . f33| . It provides 
for a nontrivial antiquark content of the nucleon at the 
starting scale of DGLAP evolution, in agreement with 
the results of global QCD fits of DIS data; see Ref. [H 
for a recent update. In particular, it quantitatively re- 
produces the flavor asymmetry of the unpolarized sea 
quarks, fi(x) — ff{x) < 0, observed in DIS [35[ and 
Drell-Yan pair production [36[ . It predicts a large flavor 
asymmetry also in the polarized sea, g±(x) — gf{x) > 0; 
there are hints of an asymmetry of this sign in a recent 
global QCD fit including semi-inclusive data [37} ; further 
clarification is ex pec ted from production in polarized 
pp collisions j38l l39l|. These nonsinglet sea quark distri- 
butions do not mix with gluons under DGLAP evolution 
and represent clear signals of the nonperturbative QCD 
vacuum structure encoded in the model. 

The matching of the model parton distributions with 
QCD quark and gluon densities is performed at the chi- 



ral symmetry-breaking scale, pr 



which represents 



the UV cutoff of the effective dynamics of constituent 
quarks. Thanks to the field-theoretical formulation of 
the dynamics (completeness of states, local interactions) 
and the relativistic covariance of the mean-field approxi- 
mation the chiral quark-soliton model conserves the over- 
all light-cone momentum, so that the constituent quarks 
and antiquarks carry the entire light-cone momentum of 
the nucleon. This provides a solid basis for matching 
the effective degrees of freedom with the quarks and glu- 
ons of QCD. Physically, the effective degrees of freedom 
are composites of the QCD quark and gluon fields, and 
the model parton distributions should be "resolved" into 



their QCD content at the scale \i 



This process is 



not governed by intrinsic properties of the effective chi- 
ral dynamics but requires detailed knowledge of its em- 
bedding in QCD, which is poorly understood at present. 
In the simplest approximation one assumes that the con- 
stituent quarks and antiquarks remain pointlike up to the 
chiral symmetry-breaking scale and matches them with 
the quarks and antiquarks of QCD; the gluon density is 
zero in this approximation [32j . This approximation was 
adopted in most calculations of partonic structure in the 
chiral quark-soliton model so far. Its accuracy may be 
judged from the fact that in leading order fits to the DIS 
data 01 at w 0.3 GcV 2 about 30% of the nucleon's 
momentum is carried by gluons. This shows that the res- 
olution effect is moderately strong in the singlet sector; 
a substantially weaker effect is expected for nonsinglets. 
More accurate matching would be possible either with 
a microscopic derivation of the effective chiral dynamics 
from QCD (such as the instanton vacuum model (24[) or 
with detailed phenomenological modeling based on em- 
pirical parton densities. 

In this article we explore the role of dynamical chi- 



ral symmetry breaking in the intrinsic transverse mo- 
mentum distribution of partons, using the chiral quark- 
soliton model as an effective description of the dynamics 
below the chiral symmetry-breaking scale. Our study is 
comprehensive and aims to address all relevant aspects of 
the problem: the definition of the pt distributions within 
the effective model, their practical evaluation and numer- 
ical study, the implementation of the UV cutoff and the 
matching with QCD, and the implications for DIS ex- 
periments. We show that the effective dynamics suggests 
a natural definition of the intrinsic transverse momen- 
tum distributions, as the momentum densities of massive 
quarks and antiquarks in the fast-moving nucleon. We 
calculate the transverse momentum distribution of va- 
lence and sea quarks in the model in leading order of the 
1/N C expansion and study their properties. Our investi- 
gation leads to several interesting new insights. 

First, we find that valence and sea quarks have very dif- 
ferent intrinsic transverse momentum distributions. The 
distribution of valence quarks (quarks minus antiquarks) 
has a range of the order of the inverse nucleon size, 
p\ ~ R~ 2 and an approximately Gaussian shape. The 
distribution of sea quarks (antiquarks), in contrast, ex- 
hibits a power-like tail oc p^ 2 that extends up to the 
chiral-symmetry breaking scale. Its coefficient is deter- 
mined by low-energy chiral dynamics and quasi model- 
independent. Such behavior is found in the flavor-singlet 
unpolarized sea quark distribution, where it was first ob- 
served in the numerical study of Ref. [4(| , and the flavor- 
nonsinglet polarized sea quark distribution, which are the 
leading combinations in the 1/N C expansion. The qual- 
itative difference between valence and sea quark trans- 
verse momenta represents the imprint of dynamical chi- 
ral symmetry breaking on the nucleon's partonic struc- 
ture and has numerous potential implications for hard 
scattering processes. 

Second, we show that, under rather general conditions, 
the sea quark transverse momentum distributions do not 
depend on the details of the UV cutoff of the effective 
chiral model. While the chiral symmetry-breaking scale 
represents the generic UV cutoff of the effective chiral 
dynamics, the manner in which it is implemented in the 
model is not dictated by chiral symmetry but must be 
constrained by other physical considerations. Imposing 
minimal physical conditions on the regularization scheme 
(charge conservation, longitudinal momentum conserva- 
tion, analyticity) we find that the sea quark transverse 
momentum distributions are independent of the regular- 
ization scheme up to momenta p\ ~ 1 GeV 2 and rep- 
resent stable predictions of the model. Since the reg- 
ularization conserves the overall light-cone momentum, 
the constituent quark and antiquark distributions in the 
model carry the entire light-cone momentum of the nu- 
cleon and can consistently be matched with QCD quarks 
and gluons at the scale p? ~ p~ 2 . 

Third, we explore the role of dynamical chiral sym- 
metry breaking at a more microscopic level, in terms of 
the light-cone wave function of the nucleon in the chiral 
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quark-soliton model [4lM43| . The large- N c limit allows 
us to discuss the nucleoli's partonic structure in terms 
of the traditional nuclear physics concepts of the mean 
field and short-range correlations. We show that the sea 
quarks in the nucleon's light-cone wave functions can ex- 
ist in correlated pairs with a transverse size of the order 
of the chiral-symmetry-breaking scale p, much smaller 
than the nucleon size R, which reflects their origin from 
dynamical chiral symmetry breaking. The pairs come 
in scalar-isoscalar (S) and pseudoscalar-isovector (II) 
quantum numbers and have a distinctive spin structure; 
at large transverse momenta p\ ~ p~ 2 their internal wave 
functions become identical, reflecting the "restoration of 
chiral symmetry" at the cutoff scale. These short-range 
correlations represent another imprint of chiral symme- 
try breaking on the nucleon's partonic structure. They 
provide a natural microscopic explanation of the high- 
Pt tails found in the sea quark transverse momentum 
distributions and point to an interesting analogy with 
short-range NN correlations in nuclei; see Refs. |44l44q | 
for a review. Most importantly, it may be possible to ob- 
serve these nonperturbative parton correlations directly 
in measurements of particle correlations between the cur- 
rent and target fragmentation regions in deep inelastic ep 
scattering or multiparton processes in pp scattering. 

Quantifying the experimental implications of our re- 
sults is a complex task, which for the most part we leave 
to a separate study. Additional information on QCD 
final-state interactions and the fragmentation process is 
needed to relate the intrinsic px distribution of partons 
to the observed transverse momentum distributions of 
hadrons emerging from hard scattering processes. Nev- 
ertheless, some simple conclusions can be drawn already 
at the present stage, without detailed modeling. Semi- 
inclusive DIS with single identified hadrons is widely 
used to measure the flavor decomposition of the nu- 
cleon parton densities. We show that the usual proce- 
dure of combining ir + and ir~ multiplicities to isolate 
the valence quark density has to be modified if the in- 
trinsic pt distributions of quarks and antiquarks in the 
nucleon are not the same and the experiment does not 
cover the full transverse momentum range of the pro- 
duced hadrons. A direct test of the nonperturbative 
parton short-range correlations predicted here could be 
performed through measurements of particle correlations 
between the current and target fragmentation regions in 
deep-inelastic ep scattering. We show that there is a 
kinematic window at moderate j*N center-of-mass en- 
ergies W 2 ~ few x 10 GeV 2 in which the two fragment- 
ing partons could be cleanly separated while pcrturbativc 
QCD radiation does not yet destroy the nonperturbative 
correlations. We also comment on the role of nonper- 
turbative correlations in multiparton processes in high- 
energy pp scattering. Finally, the nonperturbative parton 
correlations predicted here may play an important role in 
exclusive meson production at energies of W ~ few GeV. 

The plan of this paper is as follows. In Scc.[II]we sum- 
marize the model of the effective dynamics below the 



chiral symmetry-breaking scale and the resulting mean- 
field description of the nucleon in the large-7V c limit. In 
Sec. IIIII we present the definition of the transverse mo- 
mentum distributions in the model as momentum den- 
sities of constituent quarks and antiquarks in the fast- 
moving nucleon and discuss their basic properties. We 
evaluate the expressions in terms of quark single-particle 
wave functions and develop their interpretation in the 
nucleon rest frame. We also discuss the coordinate- 
space correlation function associated with the transverse 
momentum distribution in our model, and the positiv- 
ity conditions and inequalities for the polarized distri- 
butions. In Sec. IIVI we study the transverse momen- 
tum distributions of valence quarks (quarks minus anti- 
quarks) . We calculate the flavor-singlet unpolarized and 
flavor-nonsinglct polarized valence quark distributions, 

ff +d ~ u ' 3 {x,p T ) and gX- d - %+a (x,p T ), which appear in 
leading order of the 1/N C expansion, and study the aver- 
age transverse momentum (p T ). 

In Sec. [V] we give an in-depth treatment of the sea 
quark transverse momentum distributions in our ap- 
proach. We evaluate them using the gradient expansion 
of the quark Green function, an approximation which 
allows us to analytically study the behavior of the dis- 
tributions at large transverse momenta. The gradient 
expansion is formulated in terms of light-cone variables, 
which allows for a simple physical interpretation in terms 
of quark-antiquark pair production by the classical chi- 
ral field of the nucleon. We analytically exhibit the 
power-like l/p T tail of the flavor-singlet sea quark dis- 
tribution f™ +d (x,pT) and discuss its significance. We 
then describe the physical conditions on the UV cutoff 
and present two regularization schemes that meet them 
(Pauli-Villars subtraction, and an invariant-mass cut- 
off). We evaluate the distributions numerically and ver- 
ify that they are independent of the form of the UV cut- 
off. We also compute the coordinate-space correlation 
function in the model; we show that at large distances 
it decays exponentially and is completely governed by 
low-energy dynamics. Finally, we also compare the sea 
quark with the valence quark distributions at the numer- 
ical level and confirm their qualitative difference. We 
also compute the flavor-nonsinglct polarized distribution 
gi~ d {x,pT) and show that it exhibits a similar power-like 
tail at large px as the flavor-singlet unpolarized one. 

In Sec. I VII we discuss the nucleon's light-cone wave 
function at large transverse momenta. We show that it 
is dominated by configurations in which a single quark- 
antiquark pair has momenta of the order of the chiral 
symmetry-breaking scale. We compute the internal wave 
functions of S- and II-type pairs, study their spin struc- 
ture, and demonstrate that chiral symmetry is effectively 
restored at large px- We then prove that the high-py 
tails in the distribution of sea quarks, found previously 
by gradient expansion of the one-body densities, is ex- 
actly reproduced by the momentum density (wave func- 
tion overlap) of such correlated pairs. In Sec. IVIII we 
summarize the model results for transverse momentum 
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distributions and correlations and list problems meriting 
further study. We discuss the matching of the model dis- 
tributions with QCD using information about empirical 
p<r-integrated parton densities at a low scale, and discuss 
pertinent open questions. Lastly, we develop a qualita- 
tive physical picture how nonperturbative parton corre- 
lations emerge from QCD and discuss its implications. 
In Sec. I VII J we outline the implications of our results for 
hard scattering processes. We discuss the consequences 
of different p? distributions of valence and sea quarks 
for quark flavor separation in semi-inclusive DIS. We 
also investigate the possibility of probing parton corre- 
lations by measurements of particle correlations between 
the current and target fragmentation region. Lastly, we 
comment on the potential role of nonperturbative parton 
correlations in multiparton processes in pp collisions and 
exclusive meson production in ep scattering. 

A numerical study of transverse momentum-depen- 
dent quark distributions in the chiral quark-soliton 
model was reported in Ref. |4(J and found a larger aver- 
age pt of sea compared to valence quarks. Here we repro- 
duce and explain this surprising result using analytic ap- 
proximations, and show that the presence of a power-like 
tail of the px distribution follows directly from the short- 
range nature of dynamical chiral symmetry breaking in 
QCD as encoded in the effective model. We also extend 
our study to the quark helicity distribution, where the 
tail appears in the flavor nonsinglet sector and is less af- 
fected by perturbative QCD radiation in hard processes. 

Transverse momentum distributions of quarks were ex- 
tensively studied in diquark spectator models of the nu- 
cleo n |47| , the bag model [H, |49[ , light -front quark mod- 
els [5 Of. and a covariant parton model HH- While in- 
corporating some aspects of relativistic kinematics, all 
of these models describe the nucleon as system with a 
fixed number of particles, ignoring the essential many- 
body nature of the parton picture. Our approach here 
is field-theoretical and describes the nucleon as a super- 
position of configurations with different numbers of par- 
ticles, which allows us to uncover the effect of the QCD 
vacuum on the nucleon's partonic structure. Also, with 
the underlying two dynamical scales, cf. Eq. (|1.5p . our 
approach provides a parametric framework for defining 
the transverse momentum distributions at a low scale in 
terms of effective degrees of freedom, which effectively 
include also the original gauge fields of QCD. 



II. EFFECTIVE DYNAMICS AND NUCLEON 
STRUCTURE 

A. Chiral constituent quarks 

We begin by summarizing the essential elements of the 
two-scale model of low-energy dynamics and the result- 
ing description of the nucleon as a chiral soliton, follow- 
ing the lines of Refs. [3(3, IH, HH]. The effective dynamics 
resulting from the spontaneous breaking of chiral sym- 



metry can be approximated as a field theory of massive 
quarks coupled to a Goldstone boson (pion) field in a chi- 
rally invariant manner. It is described by the Lagrangian 
density 

L cS = ^(x)[td-MU^{x)]iP(x), (2.1) 

where ^p{x) is the quark field (the sum over light quark 
flavors u and d is implied), 

d = rd a , (2.2) 

and M is the dynamical quark mass. The pion field is 
contained in the variable 

U^(x) = exp[i l5 T a TT a (x)/F v ] 

= L(l +l5 )U(x) + ±(l~ l5 )UHx), (2.3) 

where 

U(x) = exp[ir a 7r a (x-)/F 7r ] (2.4) 

is the usual unitary matrix field. Here T a {a — 1,2,3) 
denote the isospin Pauli matrices, and F„ = 93 MeV is 
the pion decay constant. Typical values of the dynam- 
ical quark mass obtained from phenomcnological con- 
siderations are M ~ 0.35-0.4 GeV. In conventional 
terms, the strength of the pion-quark coupling is given 
by M/Fk sa 4, as can be seen be expanding the exponen- 
tial in Eq. <\2A\ in powers of the pion field. The effective 
theory defined by Eq. (|2.1[) is thus strongly coupled and 
has to be solved using nonperturbative methods based 
on the 1/N C expansion (semiclassical or saddle-point ap- 
proximation) [52j . 

The effective dynamics described by Eq. (|2.1[) applies 
to quarks with virtualities below the chiral-symmetry- 
breaking scale, which acts as an UV cutoff for the model. 
In practical calculations the cutoff is implemented by ap- 
plying a regularization scheme, and the actual value of 
the cutoff parameter depends on the scheme. In the fol- 
lowing we denote the generic cutoff parameter by A 2 (not 
to be confused with the QCD scale parameter Aq CD ). A 
crucial point is that the dynamical quark mass is assumed 
to be parametrically small compared to the cutoff, 

M 2 < A 2 , (2.5) 

implying that the effective model is applicable in a para- 
metrically wide range of quark momenta where the mas- 
sive quarks behave approximately as pointlike particles. 
While the numerical accuracy of this approximation is 
limited (M 2 /A 2 ~ 0.3 with a typical virtuality cutoff; 
see below), it serves as an ordering principle for the cal- 
culation of physical quantities and provides a clear math- 
ematical justification for the constituent quark picture. 
The choice of regularization scheme for implementing the 
UV cutoff involves physical judgment and is usually mo- 
tivated by the desire to preserve fundamental properties 
such as analyticity or current conservation. For the pt~ 
integrated parton densities this question was studied in 
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M 



-M 



71 classical field 
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Dirac 
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FIG. 1. (Color online) Chiral quark-soliton model of the nu- 
cleon. (a) Classical chiral field in the nucleon rest frame. The 
isospin direction at a given point is defined by the radius vec- 
tor ( "hedgehog" ) . (b) Spectrum of quark single-particle levels 
in the classical chiral field, Eq. (|2.7[) . It includes a discrete 
bound-state level and the distorted negative and positive- 
energy Dirac continua. In the nucleon state (baryon number 
+1) the bound-state and negative-energy continuum levels 
are occupied by N c quarks each. 



Refs. [32l l33f ; for the transverse momentum distributions 
of interest here it will be discussed in detail in Sec. IV Fl 
below. 

The effective model Eq. (|2.1j) can be motivated by gen- 
eral considerations about constituent quarks as effective 
degrees of freedom, see Ref. [24[ for a review. In par- 
ticular, when integrating over the quark fields and ex- 



panding the fcrmion determinant in gradients of the pion 
field, one obtains the chiral Lagrangian of the pion field, 
with definite predictions for the coupling constants in 
terms of the two parameters, M and A. In this sense, 
the effective theory "interpolates" between the chiral La- 
grangian and a theory of free quarks at large virtuali- 
tics. Eq. (|2.ip has also been derived from the instan- 
ton model of the QCD vacuum, by large-7V c bosoniza- 
tion of the instanton-induced 'tHooft interaction between 
quarks in the chirally broken phase [H, [54| • In this con- 
text the UV cutoff is given by the inverse average instan- 
ton size in the vacuum, p~ 2 « 0.4 GeV 2 , and restricts 
the quarks' Euclidean momenta (or virtualities) ; the dy- 
namical quark mass is obtained at M ~ 0.35 GeV. The 
parametric smallness of the quark mass compared to the 
cutoff, Eq. (|2.5p , follows directly from the "diluteness" of 
the instanton medium describing the QCD vacuum. 



B. Nucleon as chiral soliton 

The effective dynamics described by Eq. (|2.ip refers to 
the large-iV c limit of QCD, and information on hadron 
structure is extracted by calculating hadronic correlation 
functions in a systematic I /N c expansion. The nucleon is 
obtained as an extended system in which quarks move in 
the background of a self-consistent classical chiral field 
(relativistic mean-field approximation, or chiral quark- 



soliton model) [3Cj . In the nucleon rest frame the chiral 
field is of "hedgehog" form (see Fig.QJ,), 



U cl (x) 



exp 



-P(r) 



(r = \x\), (2.6) 



and the profile function P(r) satisfies -P(O) 
P(r) — > for r -> 



— 7r and 

55| . The quarks move in single- 



particle orbits, whose wave functions are the eigenfunc- 
tions of the Dirac Hamiltonian in the background chiral 
field, i.e., the solutions of the Dirac equation 

{-i-f-fdi+M-fUS) *„(sb) = EA(x). (2.7) 

The energy spectrum includes a discrete bound-state 
level and the positive and negative Dirac continua, dis- 
torted by the chiral field (see Fig. []J>). Each level is 
occupied with N c quarks. The static energy of the sys- 
tem is the sum of energies of the discrete level and the 
energy stored in the negative-energy continuum, 



E\U ci 



N r E\ 



lev 



N, 



Y^(E*-EW), (2.8) 



E„<0 



where En denote the energy levels of the vacuum Hamil- 
tonian with U c i = 1. The profile function P(r) is deter- 
mined by minimizing the static energy Eq. (|2.8p . In Ap- 
pendix|X]we show the profile obtained by numerical min- 
imization (5(| with a Pauli-Villars cutoff of the Dirac sea 
contribution to the energy (described in Sec.|V| and give 
a simple analytic paramctrization for use in our numeri- 
cal estimates below. The nucleon mass, which is 0(N C ) 
within the 1/N C expansion, is in leading order given by 
the minimum value of the classical energy, 



M N = E[U cl \ r 



O(N c ). 



(2.9) 



Finally, nucleon states of definite spin/isospin and linear 
momentum are constructed by quantizing the (iso-) ro- 
tational and translational zero modes of the soliton, see 
Ref. [3(| for details. This leads to the appearance of the 
N and A as rotational states of the soliton and explains 
their mass splitting, A/a — Mjsr = 0(l/N c ), as the differ- 
ence of the rotational energies. 

The description of the nucleon obtained in this ap- 
proach is fully field-theoretical and does not involve any 
ingredients "extraneous" to the dynamics encoded in 
Eq. (|2.f p . No further approximations besides the l/N c 
expansion are made in solving the dynamical problem. 
The description is also fully relativistic; it appears nonco- 
variant only because the nucleon in the large- N c limit is 
heavy, but relativistic corrections appear systematically 
as part of the 1/N C expansion. 

Matrix elements of quark one-body operators between 
nucleon states can be calculated in a systematic l/N c 
expansion and are generally expressed as sums of matrix 
elements between quark sing le-particle states in the clas- 
sical chiral field; see Ref. |31j for a review. The projection 
on nucleon spin/isospin and momentum states is done by 
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integrating over the rotational/isorotational and transla- 
tional zero modes of the classical field with appropriate 
collective wave functions. The sums over single-particle 
states can then be evaluated numerically by constructing 
the eigenfunctions and eigenvalues (both of the discrete 
bound-state level and the continuous spectrum) through 
numerical diagonalization of the Hamiltonian of Eq. (|2.7[) 
in a spherical box of finite size 57 1 . 

Alternatively, one may express matrix elements of 
quark one-body operators through the quark Green func- 
tion in the chiral background field, using the formalism 
of second quantization. In an arbitrary (generally time- 
dependent) classical chiral field, the Feynman Green 
function is defined as the solution of the inhomogeneous 
Dirac equation with a delta function source, 



[id- MlP*(x)]G F (x,y) = 6 {4) {x-y), 



(2.10) 



with causal boundary conditions, corresponding to the 
advanced solution for x° < y° and the retarded one for 
x° > y°. In second quantization this function coincides 
with the expectation value of the time-ordered product 
of quantized field operators, 

iG F (x,y) = (N\Ti/>(x)$(y)\N), (2.11) 

where \N) denotes the ground state of the fermionic sys- 
tem in the background of the (generally time-dependent) 
classical chiral field, with the occupation of the single- 
particle levels as indicated in Fig. [TJd. In the static chi- 
ral field in the nucleon rest frame, Eq. (|2.6p . the time- 
ordered product becomes 

iG F (x,y) 

= e(x°-y°) £ e- iE ^°-v ^ n (x)*i(y)j° 



n non— occ 



e(y°-x°) Y: e- iE ^°-y°^ n (x)t>l (y) 7 ° (2.12) 



2tt 



UJ - E n + iO<T n 



-1 occupied levels, 
T non-occupied levels. 



(2.14) 



Equations (|2.13p and (|2.13p can be used to convert traces 
of the Green function into sums over single-particle levels 
and vice versa. There are many practical advantages of 
working with the Green function [32|, |33| . It has simple 
transformation properties under Lorentz boosts, which 
follow directly from Eq. (|2.10[) (see Sec. IIIip . Its ana- 
lytic properties in energy allow one to derive sum rules 
and interchange sums over occupied and non-occupied 
states. Finally, the Green function can be evaluated ap- 
proximately by expanding it in derivatives of the classical 
chiral field (gradient expansion) . In this way one can ob- 
tain analytic expressions for the leading dependence of 
nucleon matrix elements in the limit of large UV cut- 
off, A — > oo, or other limiting cases, such as the parton 
distributions at large transverse momenta (see Sec. |VJ) . 



III. INTRINSIC TRANSVERSE MOMENTUM 
DISTRIBUTIONS 

A. Definition and parametric domain 

The basic framework for calculating parton densities 
in the chiral quark-soliton model of the nucleon was de- 
veloped in Refs. [32|,[33j]. Generally, in the large- N c limit 
one is interested in the parton densities at momentum 
fractions of the order 



x = OiN- 1 ), 



(3.1) 



corresponding to average (non-exceptional) configura- 
tions in the wave function of the fast-moving system [32| ■ 
The parton densities can be calculated starting from their 
parton model definition as number densities of particles 
in a nucleon moving with momentum P — > oo (33[, or 
from their representation as matrix elements of quark 
light-ray operators in the nucleon rest frame 32 1. The 
equivalence of the two formulations was demonstrated in 
Ref. [HI and is due to the relativistic and field-theoretical 
character of this description of the nucleon. 

The relativistic mean-field picture of the nucleon im- 
plies a natural definition of the intrinsic transverse mo- 
mentum distribution of partons at a low scale. Extending 
the approach of Ref. [33J, we define the intrinsic trans- 
verse momentum distribution of partons in this model as 
the pT~dependent number densities of quarks and anti- 
quarks in the fast-moving nucleon state: 



P 



(2tt) 3 
P 



(N v \ EMPMP) \Nv), (3.2) 



(N v \ EXAp)ba.(p) \N V ), (3.3) 



where pr = \pr \ 



p = (p T , xP), 
vMpf 



P = 



and the limit 



l 



is understood here and in the following. Here a aa 



(3.4) 
(3.5) 

(3.6) 

,t 



and b„ 



are the quark and antiquark annihilation 



and creation operators corresponding to the massive 
quark fields of the effective model Eq. (|2.1j) . They an- 
nihilate/create quarks in plane-wave states with three- 
momentum p and energy p° = -\/|p| 2 + M 2 , with a = 
m, d denoting the quark flavor and a the spin projection 
on the 3-axis. The notation (N v \ . . . \N V ) indicates the 
average in the ground state of the fast-moving many- 
body system, with (N V \N V ) = 1. Equations (|3.2I) and 
(|3.3[) define the unpolarizcd distributions (summed over 
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quark/antiquark spin); the corresponding expressions for 
the longitudinally polarized distributions are 



P 



(2tt) 3 
P 



(JV„| ^al a {p)a aa {p) \N V ), (3.7) 



9K^Pt) = ( N v\ &AMba*(p) \N V ), (3.8) 



where 



A a = (a = +1/2)- (a = -1/2) (3.9) 



denotes the spin difference. 

An obvious property of the transverse momentum dis- 
tributions defined by Eqs. (|3.2p and (|3.3[) is that their 
integral over px reproduce the quark/antiquark densities 
in the model, 



d 2 PT f{ l ' a (x,p T ) rc g = fi' a {x), ( 3 - 10 ) 



d Pt gi' a {x,p T ) rcg 



a. a { \ 

9i W- 



(3.11) 



The integral over transverse momenta would be logarith- 
mically divergent at large values and is rendered finite by 
the UV cutoff of the model. In our interpretation the cut- 
off is part of the definition of the transverse momentum 
distribution in the model, as indicated by the label "reg" 
in Eqs. (|3.10[) and p. lip , not as an operation applied 
only at the level of the px integral. This interpretation is 
more restrictive than the one developed in Ref. [H, (HI , 
where only the integrated parton densities were regular- 
ized, and places stronger demands on the regularization 
scheme. A detailed discussion of the implementation of 
the UV cutoff is given in Sec. IV Fl below. 

In the sense of the 1 /N c expansion the transverse mo- 
mentum we consider are of the order 



Pt 



= 0{N° C ). 



(3.12) 



The 1/N C expansion of the transverse momentum dis- 
tributions is completely analogous to that of the px~ 
integrated distributions |32j. In the following we con- 
sider the pt distributions of those spin-flavor combina- 
tions of parton densities which appear in leading order 
of the 1 /N c expansion, namely the flavor-singlet unpo- 
larizcd distributions, 

n +d {x,p T \ f! +d (x, PT ), (3.13) 

and the flavor-nonsinglet polarized distributions, 

, 9 r rf (x,p T ), 9t\x,p T ). (3.14) 

Equations (|3.2p - (|3.9p define the transverse momen- 
tum distributions of massive constituent quarks and an- 
tiquarks — effective degrees of freedom which are to be 
matched with QCD quarks, antiquarks and gluons at the 
chiral symmetry-breaking scale. The formulas should be 



regarded as preliminary definitions, whose physical sig- 
nificance and UV regularization will be elaborated in the 
following. The UV cutoff affects the distribution of va- 
lence and sea quarks very differently and will be discussed 
separately for the two cases (Secs. lTVl andlV |) . We shall see 
that under rather general conditions on the regulariza- 
tion scheme (charge conservation, longitudinal momen- 
tum conservation, analyticity) the transverse momentum 
distributions in the model are not sensitive to the details 
of the regularization scheme up to momenta p\ ~ 10 M 2 
and represent robust predictions of the model. This pro- 
vides a firm basis for the matching of the model px dis- 
tributions with QCD fSec. IVIlf . In the remainder of this 
section we discuss the formal properties of the distribu- 
tions defined by Eqs. (|3.2p - (|3.9p without reference to the 
UV cutoff. 



B. Evaluation in single— particle states 

The evaluation of the model transverse momentum dis- 
tributions defined by Eqs. (|3.2p and (|3.3p proceeds along 
the same lines as those of the pT _ integrated parton den- 
sities [33l | . In the following we summarize the main steps, 
emphasizing those aspects that require attention in the 
case of transverse momentum dependence. 

The one-body momentum density can conveniently be 
calculated in terms of the Feynman Green function of the 
system, using the formalism of second quantization. The 
quark field operator is defined as 



tp(x) 



d 3 p 



(27r)V2p° , 
b+(p)v(p,v) e** 



^2 [a<y(p) u{p,a) e" 



(3.15) 



where u(p,a) and v(p,a) are the spinor wave function 
of the free quarks and antiquarks, normalized as uu 



2M, with ■ 



W7 and v = v^j°. The quark and 



antiquark number operators can be expressed as equal- 
time products of the field operators, 



d 3 Xl d 3 x 2 e-M*!-^) ^( x ^t) il>i{x u t) 



2_, 0+ (p)a CT1 (p) — — + (3.16) 



2 Po 



d 3 xi I d A x 2 e 



— ip(xx —x 2 



E 

CT1CT2 



+ (3.17) 



where we exhibit the bispinor indices for clarity. The el- 
lipsis here represent terms corresponding to the creation 
or annihilation of particles moving in the "opposite" di- 
rection (negative momentum in the 3-direction), which 
disappear when evaluated in the many-body state with 
momentum P — > 00. It is important to note that this 
happens irrespectively of whether one integrates over the 
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transverse momentum px or keeps it fixed, and that the 
transverse momentum dependence of the structures in 
Eqs. p.l6j) and (|3.17j) surviving in the P — s- oo limit 
is unambiguously defined in our model (in conventional 
terminology these are the leading-twist distribution func- 
tions). 

The density operators required in Eqs. p.2[) and p.3j) 
can be obtained by applying appropriate spin projectors 
to Eqs. (|3.16[) and (|3.17[) . Using the standard expressions 
for the spin density matrix of a pure state (ci = a 2 = u) 
one finds that in the limit p 3 — > oo and for fixed pr 



u{p, a)u(p, a) 



(a = ±1/2). (3.18) 



Here a is defined as the quark/antiquark spin projection 
along the 3-axis. Thus the unpolarizcd density is ob- 
tained by contracting both sides of Eqs. (|3.16j) and (|3.17j) 
with 7°, 



trpy'VV'] 



r(p)<W(p), 



tr[ 7 VV^£&L(p)Mp) 



(3.19) 



(3.20) 



the longitudinally polarized density is obtained by con- 
tracting with 7° 75. 

The equal-time products in Eqs. (|3.16[) and (|3.17[) can 
be represented as appropriate limits of the time-ordered 
product of field operators, which is described by the 
Feynman Green function in the fast-moving nucleon, 

(iV,|TV(x 1 ,i 1 )^(x 2 ,t 2 )|7Y u ) = iG F {x x ,x 2 ). ( 3 - 21 ) 

Without loss of generality we choose the equal time mo- 
ment at t = and obtain 



/r +d (x, PT ) 1 Nc p 



^ j d 3 x x J d 3 x 2 e -<P(«i-«") 



f? +d (x,p T ) j (27 

x tv[G" F (0, xi; ±0, x 2 ) 7°]- (3.22) 

The Feynman Green function is given by the solution 
of the inhomogeneous Dirac equation Eq. (|2.10j) in the 
classical chiral field corresponding to the fast-moving nu- 
cleon, which is the static field of the rest frame boosted 
along the 3-direction with velocity v, 

U3"(x,t) = U£(x') 



x = 



vt 



VT 



(3.23) 



Here and in the following we use primed variables to de- 
note coordinates in the rest frame. Thanks to the Lorentz 
covariance of the Dirac equation Eq. (|2.10p and the source 
term on its right-hand side it is possible to express this 
Green function in terms of the solutions of the rest-frame 
equation. Indeed, the Feynman Green function in the 
fast-moving nucleon in the limits of Eq. (|3.22p can be 



obtained as the Lorentz boost of the occupied and non- 
occupied level parts of the rest-frame Green function, 



G F (0, Xl ;±0,x 2 ) = (±i) 



E 

n occ 

E 



/(■ 11011— occ 



x S(v) ^(x'^iix'^S^iv) e -^«(*i-4). (3.24) 

Here (t[ 2 ,s' 12 ) denote the space-time coordinates that 
transform into (0, #1.2) under a Lorentz boost in the 3- 
direction with velocity v, 



{ti 2 ,x 12 ) 



1.2 



'1,2 



(3.25) 



The Q n (x' 12 ) are the time-independent single-particle 
wave functions in the nucleon rest frame, given by the 
solution of the Dirac equation Eq. (|2.7p . and S(v) is 
the transformation matrix corresponding to the Lorentz 
boost, 

S{v) = cxp (^tV) = cosh I + 7°7 3 sinh | (3.26) 
(v = tanh?7). 



The formal proof that Eq. p.24[) represents the disconti- 
nuity of the Feynman Green function in the fast-moving 
mean field was given in Ref. [33| and relies essentially on 
the completeness of the set of single-particle wave func- 
tions. Substituting Eq. (|3.24j) into Eq. (|3.22|) one can 
straightforwardly evaluate the transverse momentum dis- 
tributions in terms of the rest-frame single-particle wave 
functions. In the limit v — > 1 one has 



v 



7 



7°+7 3 - 



(3.27) 
(3.28) 



Changing the integration variables in Eq. (|3.22p from 
X\ t2 to the rest frame coordinates x\ 2 and substituting 
the momentum representation of the rest-frame single- 
particle wave functions, 

*n«a) = / <Pv'l,2 ^ Pl2<2 ®n(p't,2), (3.29) 

one finds that they are effectively evaluated at rest-frame 
momenta 



Pi 



Pi 



(p T , ±xM n - E n 



(3.30) 



This assignment is obtained here from the constraints 
inherent in the fast-moving nucleon expression Eq. (|3.22p 
in the limit v 1; an alternative physical interpretation 
directly in the rest frame is provided in Sec. 1111 Cl below. 
Finally, substituting the nucleon momentum P as defined 
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by Eq. (|3.5[) . and taking the limit v — > 1, one obtains 

n +d {x,PT) 





f E 1 


N C M N 


n occ 


' (2tt) 3 < 


I E 1 

^ Ti non — occ ) 



[p 



X *«(p) t 7°7 + *nGP) 
(p T , ixMjv-^n)]. 



(3.31) 



[For ease of notation we have dropped the prime on the 
rest frame momenta; all momenta in the following refer 
to the rest frame unless specified otherwise (in Sec. IVI|).] 
Formally one still needs to average Eq. (|3.31[) over the 
(iso-) rotational zero modes and project on states with 
definite spin-isposin quantum numbers (32j ; however, in 
the case of the flavor-singlet unpolarized distributions 
here this operation is trivial and does not change the 
form of the expression. 

In the quark distribution in Eq. p.31[) the summation 
extends over all occupied quark single-particle states, in- 
cluding the discrete level and the negative-energy Dirac 
continuum (see Fig. [TJd). The antiquark distribution is 
given by the corresponding sum over non-occupied states 
in the spectrum, with the sign of the xM^ term re- 
versed in the 3-component of the rest frame momen- 
tum Eq. p.30|) . Equivalent representations of the quark 
distribution as a sum over non-occupied states, and of 
the antiquark distribution as a sum over occupied states, 
can be derived by making use of the completeness of the 
single-particle states (see below) [H[, 



f? +d (x,p T ) 
f? +d7 (x,p T ) 



E ) 

N c Mj\[ I n non— occ 



(2tt) 3 



E 



7 7 



"*»(p) 
[P = (pt, ±xM n - E n )\ . 



(3.32) 



In the case of the polarized distributions, the flavor- 
nonsinglct (isovector) combination appears in leading or- 
der of the 1/N C expansion. An analogous calculation as 
in the case of the unpolarized distribution leads to the 
expressions (cf. Ref.[33|) 



gr d (x,PT) ' 
gr 3 (x,PT) , 



=F 



n occ 



3(2tt) 



E 



n non— occ 



x $„(p) T T J 7 u 7 +75$r i (p) (3.33) 



[p = (p T , ±xM N - E n )] . 

Here t 3 denotes the Pauli matrix in quark flavor indices. 
The factor —1/3 arises from the integration over the 



(iso-) rotational zero modes and projection on the nu- 
cleon state with spin S3 = 1/2 and isospin T 3 = 1/2 \^ . 
An alternative representation, analogous to Eq. (|3.32|) , is 

tf-wn NM \ e 

1\ C 1VI]\{ I n non— occ 



±- 



gr d (x,PT) 



"3(2tt) 



E 



x $„(p) t r 3 7 7+75$„(p) (3.34) 

[p = Opt, ±xM n - E n )] . 

The transverse momentum distributions can also be 
expressed in terms of the Feynman Green function of 
quarks in the nucleon rest frame. It is convenient to 
introduce the energy-momentum representation of the 
Feynman Green function as 



G F {x,y) 



d 4 pi f d A p\ 



-ip 1 x+ip 2 y 



(2tt) 4 J (2tt) 4 
x2Tr6(p 1 -p°)S F (p 1 ;pi,p 2 ), (3.35) 

whose spectral representation in the rest frame is 

— $n(pl)$n(P2)h° 



^ Pi - E n + iOa n 



(3.36) 



where a n is defined in Eq. (|2.14j) . The sums over occupied 
and non-occupied single-particle levels in Eq. (|3.31[) can 
be then expressed as contour integrals in the complex en- 
ergy variable encircling the negative and positive-energy 
poles. These integrals can in turn be converted into inte- 
grals over the real energy axis by deforming the contours 
at infinity; see Ref. [32| for details. In this it is impor- 
tant that the UV cutoff respect the analytic properties of 
the unregularized theory and does not introduce spurious 
singularities. For the quark distribution one obtains 




= ±Im 



N C M N 



(2tt) 3 
x tr[S F (p°;p,p)~f 



00 df_ 

-no 27T 



(3.37) 



[p=(p T ,±xM N -p )] 



The antiquark distribution is given by the negative of 
the expression for the quark distribution with x replaced 
by —x. The corresponding expressions for the flavor- 
nonsinglet polarized distribution are 



gr d (^PT) 1 



9i 



- d (x,PT) J 



Tm ■ 



N 



dp 
, 2tt 



3(2tt) 3 
x tr[5 F (p°;p,p)T 3 7 + 7 5 ] 
[p= (p T ,±xM N -p )] . 



(3.38) 
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In the polarized case the antiquark distribution is given 
by the expression for the quark distribution with x — > — x 
without change of the overall sign. Equations (|3.37j) and 
p.38[) are the starting point for calculating the px dis- 
tribution by gradient expansion in the classical chiral 
field. They can also be used to derive equivalent repre- 
sentations of the quark distribution as a sum over non- 
occupied states, or of the antiquark distribution as a sum 
over occupied states, by deforming the integration con- 
tour in the p° integral. 




C. Rest frame interpretation 



The expression Eq. (|3.3ip for the quark and antiquark 
distributions permit a simple interpretation directly in 
the nucleon rest frame, which explains the form of the 
rest-frame three-momentum Eq. (|3.30p at which the 
single-particle wave functions are evaluated. It also re- 
veals the role of the large- N c limit in the derivation and 
emphasizes the analogy with the conventional mean-field 
approximation of nuclear physics. 

The parton densities defined by Eqs. (|3.2j) and (|3.3[) 
can equivalcntly be regarded as the densities of quarks 
and antiquarks carrying a fraction x of the nucleon's 
light-cone "plus" momentum P£ = P^ + P^. As such 
they may be evaluated in any reference frame. In the rest 

Mjv, and thus 



frame P^ 



and P^ 



n (3.39) 

(see Fig. In this frame the operator measuring the 
parton density counts quarks and antiquarks with plus 
momentum 



P 



+ = xM 



(3.40) 



In the mean-field approximation this operator removes a 
quark from an occupied single-particle state n with en- 
ergy E n and three-momentum p, which has both longi- 
tudinal and transverse components (in the case of the an- 
tiquark density, the operator puts a quark with the same 
momentum in an unoccupied state n). As the quark re- 
moved from the interacting system is not on mass-shell, 
the relation between its energy and three-momentum, 
and hence the relation between these variables and the 
external parameter x, are a priori not obvious and re- 
quire discussion. The assignment of the quark three- 
momentum of Eq. (|3.3ip obtained from the derivation in 
Sec. IIII 51 33] corresponds to 



FIG. 2. (Color online) Interpretation of the rest frame expres- 
sion of the transverse momentum-dependent quark density in 
the chiral quark-soliton model, Eq. (|3.31[) . The operator re- 
moves a quark with 3-momentum p from the single-particle 
level n in the nucleon wave function (in the case of the anti- 
quark density it would place a quark in an unoccupied level). 
The remnant system recoils with momentum p and energy 
loss — E n . The plus momentum loss calculated as the differ- 
ence of the plus momenta of the remnant system and the ini- 
tial nucleon state exactly compensates the plus momentum 
of the active quark, xMn- This shows that the relativistic 
mean-field approximation conserves plus momentum at the 
order 0(N®) (details see text). 



Because three-momentum is conserved in equal-time 
quantization, the three-momentum imparted to the rem- 
nant system in the rest frame is 



N 



O(N c ) (3.42) 



(sec Fig. [2]). At the same time, in the mean-field ap- 
proximation the energy difference between the remnant 
system and the initial nucleon state is determined by the 
energy of the removed quark, 



p o 



pO 

N 



-En 



(3.43) 



Note that the individual energies Pjj and PjJ- are 0(N C ), 
but their difference is 0(N®) and can be discussed quan- 
titatively at this level. An important point is also that 
the kinetic energy associated with the recoil of the rem- 
nant system is of the order 



E 



recoil 



p 



2M N 



OiN- 1 ) 



(3.44) 



and can be neglected relative to Eq. (|3.43l) . From 
Eqs. (|3.42[) and (|3.43p we can calculate the plus momen- 
tum difference between the remnant system and the ini- 
tial nucleon, 



P 



+ _ 



= E n +p A 



(3.41) 



We now show that this assignment can be justified by in- 
spection of the remnant system in the rest frame, i.e., the 
many-body state from which the quark or antiquark was 
removed. This system is on mass shell, and its kinematic 
variables can unambiguously be related to the external 
parameter x. 



P+ 

N 



-En 



V 



O(N c ). 



(3.45) 



This is precisely the negative of the plus momentum of 
the removed quark, Eq. (|3.4ip . obtained from the earlier 
derivation in Sec. IIII Bl Thus we see that this assignment 
implies the conservation of plus momentum of the inter- 
acting system in leading order of the 1/N C expansion, i.e., 
the individual changes in plus momenta of 0(N®) add up 
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to zero at that order. The same applies of course to the 
transverse momenta. We emphasize that this argument 
does not neglect the interactions in the system but only 
relies on the mean-field approximation to the nucleon 
wave function. The conservation of plus momentum rep- 
resents a nontrivial consequence of the consistency of the 
approximations made in this model. 



D. Positivity and inequalities 

The unpolarized transverse momentum distributions of 
quarks and antiquarks in the chiral quark-soliton model 
defined by Eqs. (|3.2|) and (|3.3|) arc particle number densi- 
ties in longitudinal and transverse momentum and should 
therefore be positive. It is easy to verify that their 
expressions in terms of single-particle wave functions, 
Eq. (|3.31[) . are indeed manifestly positive. Noting that 
the hermitean matrix 7°7 + /2 is a projector, 



7°7 H 



7°7 H 



(3.46) 



we can rewrite Eq. (|3.31l) as 

K +d (x,P T ) " 



27V c Mjv 



(2tt) 



7°7 + 



E 1 

n occ 

E 



n non— occ 



<J>n(p) 



[p = {Pt, ±xM n - E n )] 



(3.47) 



(3.48) 



Here the quark and antiquark distribution are expressed 
as a sum of explicitly positive terms. This property essen- 
tially relies on the completeness of quark single-particle 
states in the model. We note that the UV cutoff implicit 
in the effective chiral dynamics may in principle violate 
the positivity of the unregularized expression. We shall 
show in Sec. lVFI that the physical regularization schemes 
proposed there naturally preserve the positivity of the 
transverse momentum distributions. 

The completeness of quark single-particle states also 
implies the existence of inequalities between the polarized 
and unpolarized transverse momentum distributions. An 
obvious consequence of the probabilistic nature of the 
densities defined by Eqs. (j3~2]) and (|3~3| . and Eqs. ([3~7) 
and (|3.8p . is that 



f?' a (x,p T ) 



> 



\g^' a (x,p T )\ 



(3.49) 



for any given quark or antiquark flavor a [581 ]. Actually, 
in the large--/V c limit a stronger inequality was proved 
for the pT-integrated densities [59$ . namely fi(x) > 
\3gf(x)\. The generalization of this proof to the case 
of the px distributions is straightforward, and one has 



In particular, for the flavor combinations appearing in 
leading order of the 1 /N c expansion we obtain 

fl l+d (x, PT ) > \3g^- d (x,p T )\, (3.51) 

/? + <Wt)> \3gf- S (x,p T )\. (3.52) 

In Sees. IIVI and [V] we show that these inequalities are 
satisfied by the model distributions obtained here, in- 
cluding the effects of the UV cutoff. This again illus- 
trates the consistency of the scheme of approximations 
proposed here. The probabilistic character of the model 
Pt distributions expressed in the positivity condition and 
the inequalities also facilitates their matching with QCD 
quark and gluon distributions at the chiral symmetry- 
breaking scale; cf. the discussion in Sec. IVII B[ 



E. Coordinate— space correlation function 

The total (pr-integrated) parton densities in the 
chiral-quark soliton model can be represented as corre- 
lation functions of the massive quark fields of the model 
at light-like space-time separations [Hj]. Denoting the 
space-time separation four-vector by £, with £ 2 = 0, 
one usually chooses a frame where £ + = 0, £t = 0, and 
£~ 7^ and performs the Fourier transform of the cor- 
relation function with respect to £~. The equivalence 
of this representation of the parton density with the 
number density of quarks and antiquarks in the infinite- 
momentum frame was demonstrated in the chiral quark- 
soliton model in Ref. [33| . In a similar manner, the trans- 
verse momentum distributions in the model can now be 
represented in terms of correlation functions of the mas- 
sive quark fields at a finite transverse separation ^^0, 
corresponding to a space-like separation of the fields, 
£ 2 < (see Fig. [3|). It is straightforward to show that 
Eq. (|3.3ip can formally be represented as 



fl +d (x,PT) 

iT +< WT) 

fi +d {x,i T ) 1 



d2 & i(p T & 
(2^ 




(3.53) 



i r°° 

>=±-L / dC e ±ixP ^-/* 
x (Pjvl^Oh+VKO |i'jv>f+=o. 



(3.54) 

Here Pjy is the nucleon four-momentum, with ar- 



bitrary, Pj 



N.T 



0, and P, 



N 



IP, 



N 

jf. The sum 



fi ,a (x,pr) > \3g? a (x,pr)\. 



(3.50) 



over quark flavors (u,d) is implied in Eq. (|3.54j) . Con- 
versely, the coordinate-space correlation function can be 
obtained as the two-dimensional Fourier transform of the 
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T ' 



IV. VALENCE QUARK TRANSVERSE 
MOMENTUM 

A. Transverse momentum distribution 



FIG. 3. (Color online) Coordinate-space correlation func- 
tion f? +d (x,£ T ) and f? +d (x, fcr), Eq. representing the 
Fourier transform of the transverse momentum distributions 
in the chiral quark-soliton model. 



Pt distribution, 

f?"(x,Zr) = J d 2 PT e-*^ T f? +d (x,p T ) 

poo 

= 2tt / dprPT Jo(pt£t) fi +d (x,p T ), 
Jo 

(3.55) 

and similarly for the antiquark distribution. Here Jo is 
the cylindrical Besscl function. In particular, the value at 
£t = corresponds to the total (p^-integrated) parton 
density, 



r +d (x^ T = o) = .r+ d (x), 
r +d {x^ T = Q) = r + \ x ). 



(3.56) 



Here as before in Sec. IIII A[ it is assumed that the inte- 
gral over transverse momenta in the model converge; the 
role of the UV cutoff in the coordinate-space correlation 
function will be discussed in Sec. IV Gl 

We emphasize that we regard Eq. (|3.54[) strictly as 
a definition within the effective model and make no at- 
tempt to relate the bilinear quark operator in massive 
quark fields to any QCD operators with finite transverse 
separation. As discussed in Sec. HI the proper definition 
of such operators in QCD faces considerable challenges 
(gauge invariance and path dependence, rapidity diver- 
gences, renormalization properties) and is the subject of 
on-going research. In our dynamical model the opera- 
tor with finite transverse separation of the fields is well- 
defined, and its nucleon matrix element can be calculated 
within the l/N c expansion. 

Finally, the corresponding coordinate-space correla- 
tion functions for the flavor-nonsinglct polarized distri- 
bution arc 



~9i +d (x^T) 



1 



8tt _ 

x (P^|Vi(0)r 3 7+75V'(0 \Pn)s+=o- 



(3.57) 

Their relation to the px distributions is analogous to the 
unpolarized case, cf. Eq. (|3.53|) . 



Using the expressions derived in Sec. IIII Bl we now want 
to calculate the transverse momentum distributions of 
quarks and antiquarks in the chiral quark-soliton model 
and study their behavior. It will be convenient to discuss 
separately the "valence" (charge-nonsinglet, or quark mi- 
nus antiquark) and the "sea" (antiquark) distributions, 
as the two receive contributions from different classes of 
configurations in the nucleon wave function and conse- 
quently show different behavior, especially at large val- 
ues of pt- In this section we study the valence quark 
transverse momentum distributions; the sea quark dis- 
tributions will be considered in Sec. [V] 

The flavor-singlet unpolarized valence quark density 



u-\-d— u— d 



(x) = f? +d (x)-f? +d (x) (4.1) 



in the chiral quark-soliton model is dominated by the 
contribution of the discrete bound-state level in the 
quark single-particle spectrum. The contribution of the 
Dirac continuum to the valence quark density was calcu- 
lated in Ref. [32| and found to be < 5% of the total at 
non-exceptional values of x; its integral over x is zero, as 
the discrete level occupied by N c quarks already accounts 
for the entire baryon number of the nucleon. The trans- 
verse momentum distribution of valence quarks is there- 
fore also dominated by the bound-state level, and we can 
neglect the Dirac continuum contribution in our numer- 
ical study. The explicit form of the bound state wave 
function of the discrete level is given in Appendix [B] (see 
also Appendix B of Ref. [IH). Evaluating the discrete 
level term in the quark distribution given by Eq. (|3.31l) . 
and the antiquark distribution given by Eq. (|3.32[) , we 
obtain 



± 



N r M 



N 



Airp 2 
2p 3 



h 2 (p)+.f(p) 



V 



h(p)j(p) 



(4.2) 



p 3 = ±xM N 



Eley, P = \p\ 



Pt 



?,Y2 



Here h(p) and j(p) are the radial wave functions of the 
bound-state level in momentum representation as de- 
fined in Eq. (|B8[) . The valence quark pt distribution 
in our approximation is then given by [cf. Eq. (|4.1[) ] 



n 



u-\-d—u—d 



(x,Pt) = Kt d {x,PT)-m d (x,p T ). (4.3) 



We evaluate the pt distributions with M = 0.35 GeV, the 
self-consistent soliton profile described in Appendix [S] 
Eqs. (|A4| -(|A7 1) . E\ cv = 0.466 M as obtained from 
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FIG. 4. Transverse momentum distribution of flavor-singlet 
valence quarks, fi~lai( x 'P T )' Eq. (|4.1|) . as a function of p%, 
at x = 0.1 (solid line) and 0.5 (dotted line). A Gaussian 
distribution would correspond to a straight line in the coor- 
dinates of this plot. [Self-consistent soliton profile Eq. (| A4|) 
with M = 0.35 GeV and M N = 3.26 M.l 



Eq. (|B3j) . and M N = 3.26 M. As always, we put N c = 3 
when using large--/V c expressions such as Eq. (|4.2|) as a 
numerical approximation. The resulting px distributions 
are shown in Fig. [4] for two values of x. One sees that 
the valence distributions are concentrated at momenta 
of the order p T ~ M 2 , which corresponds to the inverse 
spatial size of the mean field; cf. Appendix lAl and Fig. [501 
The distributions are roughly of Gaussian shape. Such 
behavior is obtained in a variety of bound-state models 
with fixed particle number (see Ref. [sO] for a recent dis- 
cussion) and seems to be a generic consequence of the 
spatial localization of the bound-state wave function. 

The transverse momentum integral of the valence 
quark distribution is UV-finite and does not require a 
cutoff. Because the px distribution is numerically con- 
centrated at low values p T ~ few M 2 we do not consider 
the effect of possible UV regularizations on the valence 
quark px distribution. Indeed, there are physical reasons 
why the py-integrated valence quark density should not 
be regularized. Leaving it unregularized ensures that the 
baryon sum rule is satisfied in the model, as the baryon 
number resulting from the bound-state level contribution 
is unity and that from the Dirac continuum integrates to 
zero [32l l33j . It also ensures that the momentum sum 
rule is satisfied in the model HH, [56|, i.e., that the va- 
lence and sea quark distributions together carry the en- 
tire light-cone momentum of the nuclcon — a property 
of central importance for the interpretation of the model 
parton distributions and their matching with QCD (cf. 
the discussion in Sec. I VII B)) . 

Using the fact that \p 3 \ < p it is easily seen that the ex- 



FIG. 5. Average transverse momentum squared of the va- 
lence quark distribution, (pt)v&i( x ), Eq. ()4.5|l . as a func- 
tion of x. Solid line: Flavor-singlet unpolarized distribution 
ju+d-u-d rj asnec [ ij ne: Flavor-nonsinglet polarized distri- 
bution, g"" 



-d — u-\-d 



pression in brackets in Eq. (|4.2[) is positive, irrespectively 
of the sign of the radial wave functions. The discrete level 
contribution to the unpolarized quark and antiquark px 
distributions, Eq. (|4.2[) . thus satisfies 



pa 
J I, 



> o, 



< o, 



(4.4) 



for all x and px- For the quark distribution this in ac- 
cordance with the positivity condition Eq. (|3.48p . The 
antiquark px distribution in the model includes the posi- 
tive Dirac continuum contribution in addition to the dis- 
crete level contribution Eq. (|4.2p and is also positive (see 
Sec.|V|. This also guarantees positivity of the total (px~ 
integrated) antiquark density, as discussed in Ref. [32(. 



B. Average transverse momentum 

Because the distribution of valence quarks decreases 
rapidly at large transverse momenta, one may consider 
the average of p\ at fixed i as a measure of the overall 
width of the distribution, 



(Pr)vai(a;) 



J cPpx p 2 x K +d -*- d (x,pT) 
J d 2 p T f? +d - a - d (x,p T ) 



(4.5) 



Figure [S] (solid line) shows this quantity as a function of 
obtained within our approximation. One sees that 
the typical transverse momenta in the valence quark dis- 
tribution are of the order p T ~ M 2 , corresponding to the 
inverse spatial size of the mean field [cf. Eqs. (|A4[) (IA7I) ]. 
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and that the average varies only weakly with x over a 
wide range, as expected from phenomenological models 
[60| . The increase toward larger values of x indicates 
that the spatial size of configurations contributing to the 
parton density decreases at larger values of x, which is a 
general feature of a bound state with fixed particle num- 
ber. Note, however, that the mean-field approximation 
employed here is justified only for non-exceptional mo- 
mentum fractions x ~ l/N c and cannot be used to study 
the large-x limit of parton densities. 



C. Polarized distribution 

In a similar manner, the flavor-nonsinglet polarized 
valence quark density 



u — d 

= 9i 



( X )- 9 r d (x), 



(4.6) 



in the chiral quark-soliton model is dominated by the 
contribution of the discrete bound state level. The level 
contributions to the px distribution of polarized quarks 
defined by Eq. (|3.33|) . and to the antiquark distribution 
defined by Eq. (pT3i|) . evaluate to [cf. Eq. 



NqMn 
3(4ttp 2 ) 



h*(p) 



2(p 3 



— h{P)j{p) 
P 



(4.7) 



V 



±xM 



N 



E 



lev, P '■■ 



\P\ = ]/p't + (p 3 ) 2 



The resulting distribution is positive and similar in shape 
to the unpolarized distribution (see Fig. [6l dashed line) . 
The somewhat faster decrease of the polarized distribu- 
tion at large pr is due to the lower component of the 
bound-state level wave function [for j(p) = the shape 
would be the same as that of the unpolarized distribu- 
tion] and attests to the relativistic nature of the mean- 
field picture of the nucleon at large N c . The average 
(Pr) °f the polarized distribution, defined in analogy with 
Eq. (|4.5[) . is shown in Fig. [5] (dashed line) as a function of 
x. One sees that it is systematically smaller than that of 
the flavor-singlet unpolarized distribution, which again is 



a relativistic effect [61| . Similar behavior is found in rel- 
ativistic bound state models of the valence quark distri- 
butions with fixed particle number; see e.g. Refs. [49l [50| . 

It is interesting to note that the unpolarized and po- 
larized quark distributions resulting from the bound- 
state level, Eqs. (|4.2[) and (|4.7[) . satisfy the general largc- 
N c inequality for the transverse momentum distribution, 
Eq. (|3.51[) . Because numerically g^~^(x,pT) > 0, and 

also fifel(x,PT) > 0, we may replace the absolute values 
of the distributions by the distributions themselves when 
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FIG. 6. The transverse momentum distribution of fla- 
vor-nonsinglet polarized valence quarks, g±~ d ~ u+d (x,pT), 
Eq. (|4.6[) . as a function of p%, at x — 0.3. Also shown for com- 
parison is the flavor-singlet unpolarized valence quark distri- 
bution f^ +d - ii - d (x, P T) at the same value of x, cf. Eq. (14.11) 
and Fig. g] 



testing the inequality. If we then form the difference be- 
tween the left- and right-hand sides of Eq. p. 511) with 
the expressions Eqs. (|4.2[) and (|4.7[) . we obtain 



3Suev(w) 



2NM 



N 



(Anp 2 ) 



1 - 



(p 



3\2l 



P' 



J 2 (P), 



(4.8) 



which is manifestly positive because |p 3 | < p. The in- 
equality for the corresponding sea quark distributions, 
which arc dominated by the contribution of the Dirac 
continuum of single-particle quark states, is discussed in 
Sec.lVJl 



V. SEA QUARK TRANSVERSE MOMENTUM 

A. Gradient expansion 

We now turn to the transverse momentum distribu- 
tions of sea quarks in the chiral quark-soliton model. 
The sea quark distributions receive contributions from 
a broad range of quark single-particle states extending 
up to the cutoff scale. Our first concerns are to study 
how the distributions of sea quarks behave at large pt, 
how they are affected by the UV cutoff, and how they 
can be regularized in a way that satisfies basic physi- 
cal requirements (longitudinal momentum conservation, 
charge conservation, analyticity) and provides distribu- 
tions with a clear physical interpretation. We can then 
compute the sea quark px distributions numerically and 
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compare them to those of the valence quarks. For sim- 
plicity we consider first the distribution of unpolarized 
flavor-singlet sea quarks; the flavor-nonsinglct polarized 
sea will be treated summarily in Sec. IV J I 

A powerful analytic method for evaluating the sea 
quark densities in the chiral quark-soliton model is the 
gradient expansion, in which one expands the quark 
Green function in powers of derivatives of the classical 
chiral field (32[. Here we adapt this method to study 
the pt distributions. The leading-order gradient expan- 
sion contains the (exact) leading power behavior of the 
Pt distributions at large px , which reveals the role of dy- 
namical chiral symmetry breaking and allows us to study 
the effect of the UV cutoff . The leading-order expres- 
sion also provides us with an accurate "interpolating" 
approximation to the sea quark distributions at all val- 
ues of pt, which we use for a numerical estimate of the 
distributions in Sec. IV HI 

To derive the gradient expansion, we start from the ex- 
pression of the quark and antiquark pt distributions in 
terms of the Feynman Green function in the classical chi- 
ral field in the nucleon rest frame, Eq. f|3 .37[) . The Green 
function is defined as the solution of the inhomogeneous 
Dirac equation Eq. (|2.10|) . Following Ref. [62j |. we can 
regard it as the matrix element of a formal operator be- 
tween 4-dimensional "position eigenstates" |a;) = \x°,x), 



G F {x,y) = {x\[id-MU^]- 1 \y), 



(5.1) 



where id is the 4-dimensional "momentum operator" 
and d = 7"<9 Q , cf. Eq. (|2.2[) . Equivalcntly, the en- 
ergy-momentum representation of the Green function 
Eq. (|3.35| can be expressed as the matrix element of the 
same operator between conjugate momentum eigenstates 
\p) = \p°,p), with (x\p) = e-*P*, 



27rS(p 1 ~p 2 )S F (p° 1 ;p 1 ,p 2 ) 
= ( Pl \[id-MU^]- l \p 2 ). 



(5.2) 



The gradient expansion is now performed by formally 
expanding the inverse Dirac operator in gradients of the 
classical chiral field: 



[id- ML/ 75 ] -1 

= [D- 1 - MiidU 15 )}- 1 (id + MU-' 15 ) 
= [D + MD(idU^)D 
+ M 2 D(idU^)D(idU 15 )D + ..] 
x {id + MU-' 15 ), 



where 



d = (-a 2 - a/ 2 + ioy 



(5.3) 



(5.4) 



is the free massive scalar Feynman Green function. The 
leading-order contribution to the quark and antiquark 
Pt distributions Eq. (|3.37[) is contained in the terms of 
order M 2 collected after multiplying out the factors in 
Eq. (|5.3|) . Their matrix elements between 4-dimcnsional 



momentum eigenstates are calculated by inserting com- 
plete sets of momentum eigenstates between the opera- 
tors. The basic matrix elements are 

(fc 2 |E/ ±75 |fci) = 2n5(k° 1 - fc 2 °) U^{k 2 - for), (5.5) 
U*»{k) = i(l± 75 Ri(fe) 



where 



+ i(lT75)t/d(-fc) t , 



U cl (k) = / d A xe- lxk [U cl (x)-1] 



(5.6) 



(5.7) 



is the Fourier transform of the static classical chiral field 
in the rest frame [the Fourier transform of U* (x) is given 
by U(-k)i], and 

( P 2\D\ Pl ) = (2n)U^(p2-pi)D( Pl ), (5.8) 
1 



D(p) = 



p 2 - M 2 + iO' 
The relevant part of the Green function thus becomes 

Sf(p°;p,p) 



(5.9) 



M 



d 3 k 



3 [D(p)D(p-k)kU^(k)U-^(-k) 



(27T) 

[p=(p°,p), fc=(0,fc)], 



+ D 2 (p)D(p - k)kU^(k)(-k)U^(-k)p} (5.10) 



where k = k a ^ a as in Eq. (|2.2p . Substituting this ex- 
pression in Eq. (|3.37[) , evaluating the trace over spinor 
indices, and using the symmetry of the bilinear forms 
under k — > — k, we finally obtain 



4N C M N M 2 f d 3 k ~ ~ f 
= ± (2tt)3 J (2^ trfl 

x Im J ^- [D(p-k)D(p)k+ 

- D{p-k)D 2 {p)k 2 p+] (5.11) 

[p=(p°,PT,±xM N -p )]. 

Because of the symmetry of the combined momentum 
integrals the expression on the right-hand side is actually 
the same in both cases, and one has 



(5.12) 



in leading-order gradient expansion. In the following, 
for simplicity, we use the upper expression also for the 
antiquark distribution. 

Equation (|5.11l) expresses the sea quark transverse mo- 
mentum distribution in terms of the Fourier transform of 
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the classical chiral field in the rest frame and an explic- 
itly calculable loop integral over the free massive quark 
propagators, and are our starting point for the discussion 
of their physical properties. When integrating over the 
transverse momentum px, Eq. (|5.1ip reproduces the gra- 
dient expansion for the parton densities of Refs. {321 l33l|. 
The latter were shown to be equivalent to the gradient 
expansion of moments of local twist-2 operators in the 
effective chiral model. 

The gradient expansion Eq. (|5.11[) contains the leading 
power-like asymptotic behavior of the transverse momen- 
tum distribution at large momenta p\ 3> M 2 , as shown 
in detail in Sec. IV Dl below. This important property 
follows from chiral invariance, which dictates that coef- 
ficients with higher mass dimension in an expansion in 
1/Px necessarily involve higher derivatives of the clas- 
sical chiral field (63|. The gradient expansion therefore 
represents an ideal tool to evaluate the px distributions 
at large px and study the role of the UV regularization. 

In the context of the exact representation of the 
quark/antiquark distributions as sums over quark single- 
particle levels, Eqs. (|3.31|) and (|3.32p . the gradient ex- 
pansion approximates the contribution from the Dirac 
continuum, i.e., quark scattering states distorted by the 
classical chiral field. The contribution of the discrete 
bound-state level with its normalizable wave function is 
not contained in the expanded Green function. In the 
sea quark px distribution the Dirac continuum contri- 
bution is numerically dominant, especially at large mo- 
menta p T 3> M 2 , and can reliably be computed using 
the gradient expansion. At lower momenta p T ~ M 2 , as 
well as in the integral over px, the contribution of the 
discrete level becomes numerically relevant. An "inter- 
polating" approximation which includes the contribution 
from the bound-state level will be discussed in Sec. IV II 
In the following studies of the large-pr behavior of the 
distributions and the role of the UV cutoff we can safely 
neglect the level contribution and take the sea quark dis- 
tribution as defined by the gradient expansion. 



B. Representation in light— cone variables 

An interesting interpretation of the sea quark distribu- 
tion is obtained by expressing the result of the leading- 
order gradient expansion Eq. (|5.11|) in terms of light-cone 
variables. The integral over the quark rest-frame energy 
p° in Eq. (|5.11[) can be rewritten as 



dp f dp 3 



2- 



2- 



M N -p°) 

2w5{p° +p 3 -xM N ) {...} 



(5.13) 



quark 
loop 



x,p 7 
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field 
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(b) 



FIG. 7. (Color online) (a) Structure of the leading-order gra- 
dient expansion of the sea quark distribution Eq. (|5.16[) . The 
classical chiral field carries light-cone momentum fraction y 
and transverse momentum kr- It couples to a quark loop in 
which the the antiquark carries light-cone momentum frac- 
tion x < y and transverse momentum px- (b) Time-ordered 
interpretation of the quark loop integral. The classical chiral 
field creates a quark-antiquark pair with invariant mass s, 
Eq. (|5.33p . The UV cutoff of the chiral model can be imple- 
mented in the form of a cutoff in s, Eq. (|5.56ll . 



We introduce the light-cone fraction y of the momentum 
of the static chiral field (k° = 0) as 

y = k 3 /M N . (5.14) 

The light-cone components of the chiral field's momen- 
tum in the rest frame are then given by 

k+ = -k~ = k 3 =yM N . (5.15) 

Equation (|5.11[) can then be expressed as an integral over 
light-cone momentum fractions as 

'dy 



J d 2 k T M(y,k T ) 



x f q q{x,y;px,k T ) 



(5.16) 



Here f c \ denotes the light-cone momentum distribution 
of the classical chiral field of the soliton, defined as 

_ F 2 M 2 N y 



fd{y,k T ) = 



(2tt) 3 
[k=(k T ,yM N )} 



tr fl [[/ c i(fe)Uci(fe)T], (5.17) 



The function f q q describes the light-cone momentum dis- 
tribution of a quark-antiquark pair, 

t I i \ 2N C M 2 [dp- 

fqq{x,V,PT,kT) = ^fF 2 J ~2n ^~ ' 

x [D(p)k + - D 2 {p)k 2 p + ] (5.18) 
(p+ = xM N , k + = -k~ = yM N ). 
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Equation (|5.16[) permits a simple interpretation of the 
nucleon's sea quark distribution in gradient expansion 
(see Fig. [7^). The classical chiral with light-cone mo- 
mentum fraction y and transverse momentum kx "cre- 
ates" a quark-antiquark pair, of which either the quark 
or the antiquark is observed with light-cone momentum 
fraction x and transverse momentum px- Explicit calcu- 
lation below shows that f q q is nonzero only if x < y, as 
required by light-front kinematics. 

The gradient expansion Eq. (|5 . contains the leading 
asymptotic behavior of the transverse momentum distri- 
bution at large pr- For our subsequent discussion we 
would like to exhibit this behavior in a simple form. It 
turns out that in the region p\ M 2 Eq. (|5.16|) and 
its ingredients permit further simplification. The com- 
ponents of the 3-momentum of the classical chiral field 
in the rest frame are of the order of the inverse radius of 
the chiral field, 

\k T \^R~\ (5.19) 
fc 3 = yM N ~ R~ x . (5.20) 

The typical radius, as determined by the minimum of the 
classical energy of the soliton, is of the order R ~ M~ l 
[cf. Appendix [X] . For transverse momenta in the region 
p\ ^s> M 2 we can therefore neglect the dependence of 
the quark loop integral in Eq. (|5 . 16[) on kx and k~ . We 
suppose that the integral over y is dominated by non- 
exceptional values y > x, which will be borne out by 
explicit calculation (see below). Furthermore, we can 
neglect the second term in Eq. (|5.1ip . as it carries an 
extra factor of — fc 2 ~ R~ 2 . With these simplifications 
Eq. (|5.16p takes the form 

' dy 

f c \(y) fqg{x/y,p T ), (5.21) 



10 



where f c i(y) is the light-cone momentum distribution of 
the classical chiral field integrated over transverse mo- 
menta, 



Uiy) = J d\ T f cl (y, k T ) 

F 2 M 2 ii I' ~ ~ 

= - J kjjr J d " k T tr fl [U c} (k)U cl (k)^} (5.22) 

(fc 3 = yM N ). 

The function fqq(z,px) describes the light-cone and 
transverse momentum distribution of the quark-anti- 
quark pair created by the classical chiral field in the 
collinear approximation, fey = and k 2 = 0, 

2N r M 2 



fqq( Z ,PT) 



Im 



dp 
~2T 



D(p - k) 



(5.23) 



{2^fF 2 
x D(p)k + 

(p + = zk + , k" 

Here z denotes the fraction of the chiral field's light-cone 
momentum carried by the quark in the pair, 

z = x/y; (5.24) 
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FIG. 8. Gradient expansion of the sea quark transverse 
momentum distribution at x = 0.1. Solid line: Full ex- 
pression, Eq. (|5.16[) . Dashed line: Collinear approximation 
for p| > M 2 , Eq. (f5723]l . [Self-consistent soliton profile 
Eq. CM]), M = 0.35 GeV, M N = 3.26 AT.] 



the fraction carried by the antiquark is 



z — 1 — z. 



(5.25) 



We retain the quark mass M in the free propagators 
Eq. (|5.9p in Eq. (|5.23|) . Thus, formally, the approxima- 
tion Eq. (|5.21[) corresponds to the limit of large radius of 
the chiral field, R S> M _1 , at fixed quark mass M. We 
refer to Eq. (|5.21j) as the collinear approximation. It has 
the same structure as the DGLAP evolution equations 
describing parton splitting in perturbativc QCD in the 
collinear approximation. Note that the simplifications 
made here — integration over the chiral field's kr inde- 
pendently of pt, and neglect of k 2 — are parametrically 
justified in the domain p\ M 2 , and that Eq. (|5.21[) 
captures the exact leading behavior of the model px dis- 
tribution in this region. 

In Fig. [8] we compare the pr distributions obtained 
from the full expression Eq. (|5.16[) and the collinear ap- 
proximation Eq. (|5.23[) at a typical values x = 0.1 (the 
distributions shown here were regularized with a Pauli- 
Villars cutoff, described in Sec. IV Fl below). One sees 
that the collinear approximation accurately reproduces 
the full gradient expansion with better than 30% accu- 
racy above p\ ~ 5 M 2 . In the following investigations of 
the structure of the distributions at large px and the role 
of the UV cutoff we can therefore use the simpler collinear 
expression Eq. (|5.16[) ; for numerical estimates at finite pt 
we use shall use the full expression Eq. (|5.16[) . 
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FIG. 9. Light-cone momentum distributions of the classical 
chiral field in the large-iV c nucleon. Solid line: Flavor-singlet 
unpolarized distribution f c i(y), Eq. (|5,22p . as appears in the 
convolution integral Eq. (|5.21|) . The distribution shown here 
corresponds to the self-consistent soliton profile Eq. (|A4fl . 
with M = 0.35 GeV and M N = 3.26 M. Dashed line: Flavor- 
nonsinglet polarized distribution g c i(y), Eq. (|5.7ip . as appears 
in Eq. (|5.70[) . The plot shows 3g c \(y), which is constrained 
by the positivity condition Eq. (|5,78[) . 



D. Structure of quark— antiquark pair 

We now turn to the light-cone momentum distribution 
of the quark-antiquark pair in the gradient expansion of 
the sea quark distribution, denned by Eq. (|5.18|) . We first 
study it without the UV cutoff; the implementation of 
the cutoff will be considered in the following subsection. 
The integral over p~ can be calculated straightforwardly 
by closing the integration contour in the complex plane. 
The poles in p~ are located at 



p — k 



(pr - k T f + M 2 - iQ 



Pt 



p+-k+ 
M 2 - iO 



(5.26) 
(5.27) 



The integral is nonzero only if they lie on opposite sides 
of the real axis, which requires k + > \p + \. In terms of 
the variable 



p + /k + , 



(5.28) 



which measures the fraction of the pair's light-cone mo- 
mentum carried by the quark, cf. Eq. (|5.24[) . this condi- 
tion implies \z\ < 1. Because x > we can limit ourselves 
to 2 > 0, and thus < z < 1. We close the contour 
around the pole Eq. (|5.26j) in the upper half plane. At 
the pole, the virtuality of the active quark is 

p 2 — M 2 =t — M 2 

= --[(pr - zk T f + M 2 - zzk 2 }, (5.29) 



Momentum distribution of classical field 



We first consider the light-cone momentum distribu- 
tion of the classical chiral field, defined by Eq. (|5.22p . 
It can be evaluated using the explicit expression given 
in Eq. (|C5|) in Appendix [D] and the soliton profile 
paramctrization Eqs. (|A4[) - (|A7[) in Appendix [A"l The nu- 
merical distribution is shown in Fig. [5] (solid line). One 
sees that the light-cone momentum distribution of the 
chiral field extends over a broad range of momentum frac- 
tions y ^ 0.1 — 0.5. It is not limited to values (/ « 1, 
as would be the case for single pions emitted by valence 
quarks, and in this sense reflects the complex interactions 
in the nucleon in the mean-field approximation (cf. our 
discussion of parton correlations in Sec. IVIIBI below). 
Note that in this classical approach based on large- N c 
limit the chiral field's light-cone momentum fraction y 
is not limited to values y < 1; however, the distribu- 
tions become exponentially small in the limit y — > oo 
}32l ]. Comparison of the results obtained with the self- 
consistent soliton profile Eq. (|A4[) and the variational 
profile Eq. (|A8|) shows that the bulk distributions are 
not sensitive to the details of the profile and represent 
stable features of the model 64 1. 



where we have used that k + k = k 2 + k 2 -,. The integral 
in Eq. (|5 . 18[) then becomes 



fqq(x,y,PT,k T ) = 



2N C M 2 0(z)9(z) 
1 zk 2 



M 2 -t (M 2 -tf 



(5.30) 



We can equivalently express the quark-antiquark mo- 
mentum distribution in terms of the variables of light- 
front time-ordered perturbation theory (see Fig. [TJd). 
In this approach light-front 3-momenta (i.e., the plus 
and transverse components of the 4-momcntum) are con- 
served and intermediate particles are on mass shell, but 
the light -front energy (i.e., the minus component of the 
4-momentum) of the intermediate state is generally dif- 
ferent from that of the initial state, resulting in non- 
conservation of 4-momentum. With the (+,—,_!_) com- 
ponents of the 4-momenta of the quark and antiquark 
given by 



, p 2 T + M 2 
I zk . zk + » Pt 



P2 



(p T - k T ) 



M 2 



zk^ 



(5.31) 



hr -pr), (5-32) 
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the invariant mass of the quark-antiquark pair is 



s = (pi 



P2) 



(p T - zk T ) 



M 2 



(5.33) 



This corresponds to the invariant mass of a pair with zero 
overall transverse momentum, subjected to a transverse 
boost by fc<r- Thus, we find the following simple rela- 
tion between the active quark virtuality in the invariant 
approach and the invariant mass difference in the time- 
ordered approach: 



i-M 2 = -z(s-k 2 ). 



(5.34) 



The quark-antiquark distribution Eq. (|5.30l) then takes 
the form 



fqq(x,y,PT,k T ) 



2N C M 2 Q(z)Q(z) 



(2nfF 2 
1 



zz 
k 2 



-k 2 (s-fc 2 ) 2 
2N C M 2 Q(z)Q(z) s 
(2tt) 3 F 2 ii (s-fc 2 ) 2 ' 



(5.35) 



In the collinear approximation Eq. (|5.23[) , if we neglect 
the overall transverse momentum kx and the virtuality 
fc 2 according to Eqs. ([5TT9"]) - ([5T2"0]) . the quark-antiquark 
momentum distribution becomes 



fqqi Z ,PT) = 



2N C M 2 Q(z)@(z) 



(2tt) 3 F 2 p 2 T + M 2 



(5.36) 



This result has several interesting features. First, the 
transverse momentum distribution in the pair exhibits a 
power-like "tail" at large values, 



fqq{z,PT) 



2N C M 2 9(z)9(f) 



Pt 



{p 2 T > M 2 



{2nfF 2 

(5.37) 

Through Eq. (|5.21l) it produces a similar tail in the nu- 
cleon's sea quark transverse momentum distribution; this 
feature will be discussed in detail in Scc. lVEl bclow. Sec- 
ond, the ~ tail results in a would-be logarithmic di- 
vergence of the quark-antiquark density in the pair when 
integrated over pt up to the cutoff scale. The coefficient 
of this logarithmic divergence is the same as that of the 
logarithmic divergence of F 2 in the effective theory, 



Fz 



N C M 2 



A 



2 l0g i^ 



(A 2 > A/ 2 ). (5.38) 



As a consequence, in the limit A — > oo the quark- 
antiquark distribution in the pair is normalized such that 



A 2 



d 2 PT fqq{z,p T ) 



1 



(A 2 -> oo). (5.39) 



We shall require that this condition be satisfied also for 
finite values of the cutoff when we introduce the UV reg- 
ularization in Sec. IV Fl 



We note that in Eq. (|5.36[) the quark-antiquark mo- 
mentum distribution in the pair is actually independent 
of the momentum fraction z. This is because at this 
level of approximation there is no restriction on the con- 
stituent quarks' virtuality or invariant mass. A nontriv- 
ial z-dependence of the distribution will appear when a 
cutoff is introduced. 



E. Power— like tail at large momenta 



Through the convolution formula of Eq. (|5.2ip the 
l/p T tail in the momentum distribution of the quark- 
antiquark pair, Eq. (|5.37[) . causes a similar behavior in 
the nucleon's sea quark transverse momentum distribu- 
tion: 



1 



Cl +d (x) 
Pt 



(p|»M 2 ), (5.40) 



fi+d - _ 2N C M 2 



— /d(y)- 
v 



(5.41) 



The pt distribution of sea quarks thus has a power- 
like behavior at momenta p T ^> M 2 . This fact is of 
central importance and implies that the distribution of 
sea quarks is qualitatively different from that of valence 
quarks (see Sec. HV|) . It is a direct consequence of dynam- 
ical chiral symmetry breaking and shows the imprint of 
the QCD vacuum on the nucleon's partonic structure. 
Note that the coefficient of the tail is determined by the 
effective chiral dynamics at the scale ~ M and can be 
computed without explicit reference to the UV cutoff of 
the effective theory. In fact, the ratio M/F^ appearing in 
Eqs. (|5.37|) and (|5.41[) is just the coupling constant of the 
massive constituent quarks to the chiral field, as it fol- 
lows from expanding the effective interaction Eq. (|2.ip . 
cf. Sec. Ill Al A corresponding interpretation of the tail 
as due to quark-antiquark pairs created by the classical 
chiral field will be developed in Sec. IVII 

Paramctrically, the power-like tail in the sea quark pt 
distribution extends up to the UV cutoff scale A 2 . It 
is interesting to note that the coefficient of the tail re- 
lated to the would-be logarithmic divergence of the pt~ 
intcgrated parton density and can be recovered as the 
derivative of the pT~integrated parton density with re- 
spect to the upper limit of the pt integral, in formal 
analogy to the relation for the unintcgrated parton den- 
sity in perturbative QCD, Eqs. (|1.2p and (|1.3p . In fact, 
the coefficient Eq. (|5.4ip is nothing but the pt integrated 
parton density in gradient expansion, 



r it 
J 1 



u+d / \ 

;radW- 



(5.42) 



However, we caution that Eqs. (|5.40p - (|5.42p apply with 
logarithmic accuracy only, and that the numerical values 
of the pt distribution for finite cutoff arc strongly af- 
fected by the cutoff. In Scc. lVFl wc formulate the physi- 
cal conditions under which the UV regularization should 
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be implemented, and show that the high-momentum 
tail of the sea quark distribution — albeit in numerically 
modified form — is indeed a robust feature of the model. 



F. Implementation of ultraviolet cutoff 

To proceed further with our study of sea quark trans- 
verse momentum distributions we now have to discuss 
the implementation of the UV cutoff of the model and 
its effect on the distributions. This will allow us not only 
to evaluate the distributions quantitatively, but also to 
integrate them over pr and relate them to the total par- 
ton densities. To study the effects of the UV cutoff we 
use the gradient expansion in the collinear approxima- 
tion, Eq. (|5.2ip . which captures the leading behavior of 
the unregularized distributions at large pr and allows us 
to illustrate the essential points in analytic form. 

As explained in Sec. HI the manner in which the cutoff 
is implemented defines the effective degrees of freedom of 
the model and must be based on physical considerations 
going beyond the intrinsic structure of the effective chiral 
theory. Here we require that the regularization procedure 
satisfy the following conditions: 

I) The regularized distribution should preserve the 
symmetry of quarks and antiquarks in the pair: 



fqq(z,PT)re S = fqq 0, PT^eg- 



(5.43) 



Exchange of quark and antiquark amounts to ex- 
changing z — > z = 1 — z and pr — > — Pt', because 
the unpolarized distribution is a function only of 
Pt = \pt\ the latter change is trivial. 

II) The regularized distribution should be normalized 
such that the total number of quarks and anti- 
quarks in the pair is unity, and that they carry the 
entire longitudinal momentum of the chiral field. 
This implies that the pT~integrated distribution 



fqq( z )rc 



satisfy 



! = J d 2 p T fqq(z,Pr] 
/ dz fqq(z) Icg = 1, 

Jo 

f 1 1 

J dz Z fqq(z) Icg — —■ 



(5.44) 



(5.45) 
(5.46) 



Because of the symmetry with respect to z — > 
1 — z the two requirements are actually identical. 
Namely, 



dz Z fqq(z^) re g 



dz z /, 



qq\ z )i-cg 



1 f 1 

- / dz (Z + Z) f q q(z) rcg 
1 



dz fqq{z) 



rog- 



(5.47) 



Physically, these conditions imply that the massive 
quarks and antiquarks are the only constituents of 
the nucleon's light-front wave function up to the 
cutoff scale, and that there is no momentum "hid- 
den" in other degrees of freedom. 

Ill) The regularization should not change the large- 
distance behavior of the quark field correlator in co- 
ordinate space. This requirement will be discussed 
in detail in Sec. IV Gl and implies that the cutoff 
should not modify the analytic properties of the px 
distribution in the vicinity of the leading singularity 
in pt at complex values of the order M, which gov- 
erns the fall-off at large distances. In the collinear 



approximation this is the pole at p\ 
Eq. (15361 . 



The above represents a minimal set of physical require- 
ments based on our present understanding; they may be 
amended by further conditions if more insights into the 
nature of the effective degrees of freedom were gained in 
the future. In the following we explore to what extent 
these minimal requirements determine the sea quark pt 
distributions quantitatively. We first present two regu- 
larization schemes that meet these requirements. 

Pauli-Villars subtraction. In the Pauli-Villars (PV) 
regularization scheme [HI, [33| one subtracts from the 
original distribution of pointlike quarks with mass M a 
certain multiple of the analogous distribution of quarks 
with a regulator mass Mpy > M , 

fqq( Z ,PT)pV = fqq(z,p T \M) 
M 2 

~ T72-frt(*,PT\Mrv). (5.48) 

The coefficient is chosen such that the subtraction cancels 
the logarithmic divergence associated with the integral 
over pt- Applying this subtraction to Eq. (|5.36|) we get 

2N M 2 

fqq(z,PT)pV = ^J^ ^ Q{z)Q(z) 



(2nfF 2 



M 2 W - M 2 



(p 2 T + Mi)( P 2 T + M 2 y y 



(5.49) 



The regularized distribution drops as ~ 1 jp T at p? — > oo 
and is integrable. At the same time, we replace F 2 in 
the normalization factor by the result obtained with PV 
regularization, 



F 



N C M 2 



lO! 



Mpv 
M 2 ' 



(5.50) 



One can easily verify that with this normalization the 
conditions Eqs. (|5.45[) and (|5.46|) arc satisfied, and that 



fqq(z) 



PV 



= j d 2 p T fqq(z,PT)pV 



1. 



(5.51) 



Note also that the subtraction docs not change the 
residue of the pole at p\ = ~M 2 \ i.e., Eq. (|5.49|) has 
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the same behavior near p T — > —M 2 as the unregular- 
ized expression Eq. (|5.36[) . The numerical value of the 
regulator mass is fixed by requiring that Eq. (|5.50p re- 
produce the physical value of the pion decay constant, 
F„ = 93 MeV. One obtains 



M| v /M 2 = 2.52 (M = 0.35GeV). 



(5.52) 

Invariant mass cutoff. Another way of implementing 
the cutoff is to limit the invariant mass of the quark- 
antiquark pair in the time-ordered approach, Eq. (|5.33[) 
(see Fig. [TJd) . In the collinear approximation the invari- 
ant mass Eq. (|5.33j) is given by 



Pt 



M 2 



> 0. 



(5.53) 



The cutoff is implemented by multiplying the vertices of 
the quark-antiquark pairs in the initial and final state 
with a form factor F(s) that suppresses invariant masses 
of the order s ~ A 2 . Wc consider a monopole and an 
exponential form factor, 



F(s) 



1 



1 + s/A 2 ' 
F(s) exp = exp(-s/A 2 ). 



(5.54) 
(5.55) 



The form factor is normalized to unity at the unphysical 
point s = 0, which corresponds to p T = — M 2 ; the signif- 
icance of this choice will be explained in Sec. IV Gl The 
regularized quark/ antiquark distribution is then given by 

fqq(z,PT)inv = fqq(z,p T ) F 2 (s) 

2N C M 2 0(z)9(z) F 2 (s) 
~ (2tt) 3 F 2 p 2 T + M 2 ' 

The normalization condition now takes the form 

F 2 (s) 



(5.56) 



2N C M 2 
WfF 2 



dz I d pt 



Pt 



M 2 



1. 



(5.57) 



The value of A 2 is fixed such that Eq. (|5.57|) is satisfied 
with the physical value of F 2 . For M = 0.35 GeV this 
gives 



A 
A 



2 

mon 
2 

exp 



31.1 M 2 
44.8 M 2 



7.78 x 4Af 2 , 
11.2 x 4M 2 . 



(5.58) 
(5.59) 



The latter values are given as multiples of the minimum 
value of s in the physical region, 4M 2 . Figure fTUk shows 
the effective z-distribution obtained with an invariant 
mass cutoff (exponential form factor) for several values 
of pt] for the sake of comparison all distributions were 
normalized to unit integral over z. One sees that the 
z-distribution is rather fiat for pt = and becomes pro- 
gressively more concentrated around z = 1/2 as pt in- 
creases. 

The px~mtegrated distributions in the quark-anti- 
quark pair, Eq. (|5.44[) . obtained with the two rcgulariza- 
tion schemes, are shown in Fig. 110b . Several features are 
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FIG. 10. (a) Quark-antiquark light-cone momentum distri- 
bution f q q(z,pT) as function of z for fixed pr [exponential 
invariant mass cutoff, Eq. (|5.55p ]. The values of are in- 
dicated above/below the curves. For the sake of comparison 
the distributions at all values of p% were normalized to unit 
integral over z. (b) The p-r-integrated quark-antiquark light- 
cone momentum distribution f qq -(z), Eq. I|5.44|) . with differ- 
ent UV cutoffs. Dotted line: PV subtraction. Dashed line: 
Invariant mass cutoff by monopole form factor, Eq. (15.5411 . 
Solid line: Invariant mass cutoff by exponential form factor, 
Eq. 1(5351) . 



worth noting: (a) The light-cone momentum distribu- 
tion obtained with the PV cutoff is independent of z and 
given by fqq(z) = 1. In the context of the convolution 
integral, Eq. (|5.21[) , this turns out to be a reasonable ap- 
proximation for x that are not parametrically small (see 
below), (b) The light-cone momentum distribution ob- 
tained with the invariant mass cutoff vanishes at the end 
points z — > 0, 1. In the convolution integral Eq. (|5.21l) . 
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FIG. 11. Sea quark transverse momentum distribution 
fi~grtd( x ' P T ) calculated using the gradient expansion in 
collinear approximation, Eq. (|5.21[) , as a function of p?, at 
x — 0.1. Shown are the distributions obtained with var- 
ious UV cutoff schemes (see text). Both pr and /"g^d 
are given in units of the constituent quark mass M. Dot- 
ted line: Unregularized distribution (no UV cutoff). Solid 
line: PV subtraction. Dashed line: Invariant mass cutoff, 
exponential form. Dash-dotted line: Invariant mass cutoff, 
monopole form. [Self-consistent soliton profile Eq. (|A4|) with 
M = 0.35 GeV and M N = 3.26 M] 



dashed, and dash-dotted lines. The results show several 
notable features. First, the distributions with any UV 
cutoff differ from the one without cutoff already at low 
values p T ~ few M 2 . It shows that the cutoff plays an es- 
sential role in the numerical value of the px distribution 
already at low px- This fact is not obvious from para- 
metric considerations based on the hierarchy A 3> M, as 
the distributions at fixed pt are UV finite and thus do 
not "require" regularization. 

Second, the distributions obtained with PV subtrac- 
tion and the invariant mass cutoff (exponential and 
monopole) are very close up to values p T sa 10 M 2 . In 
this region of px they are determined by generic features 
of the cutoff as are fixed by our general conditions. This 
finding is very important, as it means that the pt distri- 
butions in this range are robust predictions of the model 
and can be discussed at the same level as other low- 
energy characteristics of the nuclcon. 

Third, at p\ > 10 M 2 the distributions obtained with 
the different cutoffs show large differences, as expected. 
One important implication of this is that the averages 
differ substantially and do not serve as reliable 
measures of the width of the bulk of the pt distribu- 
tion. These averages assign disproportionate weight to 
the high-j»T region where the distributions are not con- 
strained by our requirements on the cutoff. 



G. Quark field correlator in coordinate space 



because f c \(y)/y is finite in the limit y — > 0, the van- 
ishing of fqq{z) at z —> ensures that f(x,px) — > for 
x — > 0, i.e., the sea quark density in the nucleon vanishes 
at x — > 0. (c) With both cutoff schemes the convolution 
integral for the sea quark transverse momentum distri- 
bution Eq. (|5.21[) involves a broad distribution of quark 
momentum fractions centered around z ~ 1/2, at least 
for non-exceptional values of x. 

With all the ingredients in place, we can now calculate 
the sea quark transverse momentum distributions for a 
finite cutoff and study their dependence on the UV regu- 
larization. We use the gradient expansion in the collinear 
approximation, Eq. (|5.21[) . with the light-cone momen- 
tum distribution of the chiral field given by Eq. (|5.22p and 
evaluated with the soliton profile of Eqs. (fA3 )) -(|A7 |l [cf. 
Fig. [9]. The px distribution in the quark-antiquark pair 
is now given by the regularized expressions Eq. (|5.49[) 
or Eq. (|5.56[) . The numerical results for the px distri- 
butions at a representative value of x = 0.1 are sum- 
marized in Fig. 1111 The py-distribution without cutoff 
(or, what is the same, in the limit of infinite cutoff), 
as obtained from the unregularized quark-antiquark dis- 
tribution Eq. (|5.36[) . is shown by the short-dashed line; 
this distribution is not integrable over px and shown 
for comparison only. The distributions obtained with 
the PV cutoff and the invariant mass cutoff (exponential 
and monopole form factor) are shown by the solid, long- 



Further insight into the role of the UV cutoff in the sea 
quark transverse momentum distributions can be gained 
by studying the corresponding coordinate-space correla- 
tion functions (cf. Sec. IIII Ep . The essential points can 
again be illustrated using the gradient expansion in the 
collinear approximation, Eq. (|5.21[) . In this approxi- 
mation the coordinate-space correlation function of sea 
quarks, Eq. (|3.54[) . is given by the transverse Fourier 
transform of the convolution formula, Eq. (|5.21[) , 

/? +< Wt)= f°°- Uy)f g g(x/y^ T ), (5.60) 

j x y 

where f q q(z, £t) (z = x/y) is the transverse Fourier trans- 
form of the px distribution in the quark-antiquark pair, 

f qq (z,fr) = J d 2 PT e-*^ f qq (z^ T ) 

p OC 

= 2?r / dpTPT Jo(pt£t) fqq(z,p T ). (5.61) 
Jo 

This function describes the spatial structure of the 
quark-antiquark pairs created by the chiral field. It is 
interesting to see how it behaves at small and large dis- 
tances, and how its behavior is modified by the UV cutoff 
of the model. In particular, the behavior at large dis- 
tances sheds new light on the regularization conditions 
put forward in Sec. IV Fl 
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FIG. 12. Coordinate-space correlator /"^.^(^^t), corre- 
sponding to the gradient expansion of the sea quark distri- 
bution in collinear approximation, Eq. (|5.60|) . as a function 
of £t, at x = 0.1. Dotted line: No cutoff. Dashed line: 
PV subtraction. Solid line: Invariant mass cutoff (exponen- 
tial). [Self-consistent soliton profile Eq. (|A4| with — 0, 
M = 0.35 GeV, M N = 3.26 M.j 

Without a cutoff, i.e., with the distribution Eq. (|5.36l) . 
the Fourier transform of the pr distribution in the qq pair 
is given by 

U{zAt) = 2tt Ko(M£ t ). (5.62) 
At £t — t this function diverges logarithmically, 

f qq (z, fr) ~ 2tt log J- (Ct -> 0), (5.63) 

IVl l^T 

which reflects the logarithmic divergence of the total (pt~ 
integrated) parton density in the model without cutoff. 
With a cutoff this divergence is regularized. For example, 
with the PV subtraction the distribution becomes 

/m(*,6t)pv = 2tt [K (Mfr) - Ko(Mpv^t)} 
Mpv 

~27rlog^ «r->0). (5.64) 

A similar result is obtained with the invariant mass cut- 
off. Thus we see that at small distances the behavior of 
the correlation function is qualitatively changed by the 
cutoff, in accordance with the fact that the function at 
£t = coincides with the total parton density 

At large distances the £t correlation in the qq pair 
decays exponentially, 

~2tt [V(2M£t)] 1/2 

x exp(-Af£ T ) (£ T oo), (5.65) 

with a range determined by the constituent quark mass 
M . This appears natural, as the constituent quark mass 



represents the "mass gap" of the effective chiral model 
and controls the long-distance behavior of quark corre- 
lations. Since it reflects a low-energy property of our 
model we should require that the UV cutoff do not mod- 
ify this behavior. The exponential fall-off at large dis- 
tances is related to the singularity of the unregularized 
Pt distribution at pip = — il/ 2 , corresponding to complex 
values oipT- Our requirement therefore implies that the 
cutoff should not modify this free-field singularity. It is 
easy to see that the schemes discussed in Sec. IV Fl satisfy 
this requirement. PV subtraction leaves the residue of 
the pole at p\ = —M 2 unchanged, cf. Eq. (|5.49[) . Like- 
wise, with the invariant mass cutoff the residue remains 
unchanged because the invariant mass vanishes at the 
pole, s — > 0, and the form factors are normalized such 
that F(s = 0) = 1, cf. Eq. ([535]) . This explains the 
physical basis of regularization condition III presented in 
Sec. EH 

The quark field correlator fi +d (x,^T) m the nucleon 
obtained from Eq. (|5.60[) is presented in Fig. [12j The 
plot shows the £t dependence of the correlation function 
at a representative value of x = 0.1. The following fea- 
tures are worth noting: (a) The regularized distributions 
rapidly approach exponential behavior at £r ~ fewM , 
as implied by Eq. (|5.65[) . (b) The distributions obtained 
with the PV and invariant mass cutoffs arc overall rather 
close at all distances. This explains the approximate 
cutoff-independence of the p^-distributions at low px 
observed in Sec. rVTl (see Fig. [TTJ . We note that the 
results would be even closer if we required that the dif- 
ferent regularizations reproduce the same total parton 
density, e.g. by adjusting the cutoff values, or by per- 
forming a second PV subtraction, (c) The curvature of 
the £t distributions at £t = 0, which is proportional to 
(Pt) i is n °t effectively constrained and can thus vary con- 
siderably between different regularization schemes. This 
reaffirms our earlier conclusion regarding the model de- 
pendence of (p T ). 

In sum, we find that the combined requirements of 
matching the parton density at £t = and decay- 
ing exponentially at £t - > oo effectively constrain the 
correlator at all distances. The Fourier transform of 
these correlators results in stable, cutoff-independent 
px-distribution at low px as seen in Fig. 1111 The study of 
the coordinate-space correlation function of sea quarks 
thus reaffirms our conclusion that the physical regular- 
ization conditions presented in Sec. IV Fl result in robust 
transverse momentum distributions at low px- 

Our aim in this section was to investigate the influence 
of the UV cutoff on the coordinate-space distribution, 
which could be done in a simple way with the gradient 
expansion in collinear approximation, Eq. (|5.21[) . In this 
approximation the information about the finite spatial 
size of the classical chiral field is lost, and the exponen- 
tial decay of the correlator at large £t is due entirely to 
the intrinsic size of the quark-antiquark pair. To study 
the "true" large-distance behavior of the sea quark cor- 
relator in the nucleon one should use the full gradient 
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expansion result, Eq. (|5.16[) . in which the finite size of 
the chiral field is encoded in the k T dependence. Numer- 
ical studies show that also in the full correlation function 
the large-^T behavior is independent of the UV cutoff; 
i.e., our conclusions are general and do not depend on 
the additional simplifications made in Eq. (|5.21|) . 



H. Numerical evaluation 

Having established the behavior of the sea quark trans- 
verse momentum distribution at large momenta p\ 3> 
M 2 and the role of the UV cutoff, we now want to make 
a numerical estimate of the distributions also at lower 
momenta p\ ~ few M 2 . This will allow us to compare 
the sea quark transverse momentum distribution with 
those of the valence quarks at a quantitative level (see 
Sec. EH). 

In its representation as a sum over occupied quark 
single-particle levels, Eq. (|3.32j) . the sea quark distribu- 
tion receives contributions both from the negative-energy 
Dirac continuum and the discrete bound-state level. The 
gradient expansion approximates the contribution from 
the Dirac continuum, i.e., quark scattering states dis- 
torted by the classical chiral field, which dominates at 
momenta p\ 3> M 2 . The contribution from the discrete 
bound-state level with its normalizable wave function is 
not contained in the expanded Green function. At lower 
momenta p\ ~ fewM 2 , as well as in the integral over 
Pt, the contribution of the discrete level becomes nu- 
merically relevant. A more accurate approximation in 
this region for numerical purposes is obtained by adding 
to the leading-order gradient expansion Eq. (|5.21j) the 
contribution from the discrete bound state level given 
by Eq. |jO) . 

/? +< V,Pt) « f^Li^PT) + f?${x iPT ). (5.66) 

This approximation is known as the "interpolation for- 
mula," as it becomes exact both in the limit of large 
soliton size, where the gradient expansion is paramet- 
rically justified and the discrete level energy becomes 
negative, and in the limit of small soliton size, where 
the level contribution dominates [HI]. Numerical stud- 
ies show that Eq. (|5.66p reproduces the exact numerical 
result for the disintegrated sea quark distribution (56j 
with an accuracy of far better than 20% for x = 0.1 — 0.5 
when evaluated with the self-consistent soliton profile 
Eqs. (|A4[) - (|A7p . We therefore expect it to provide a rea- 
sonable approximation also for the px distributions at 
p\ ~ few M 2 . 

Figure [T3l shows the different contributions to the sea 
quark distribution in Eq. (|5.66|) at a representative value 
of x = 0.1 (note that the level contribution is nega- 
tive and shown with opposite sign in the figure). One 
sees that the gradient expansion contribution approxi- 
mating the Dirac continuum clearly dominates at large 
p T ; the level contribution to the total ff +d is < 20% 
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FIG. 13. Numerical approximation for the sea quark trans- 
verse momentum distribution fi +d (x,pr), Eq. (|5.66[1 ("in- 
terpolation formula"). Dashed line: Gradient expansion 
fl^fad( x iP T )' Eq. <|5.21[) . approximating the Dirac continuum 
contribution in Eq. ()3.32[) (PV regularization). Dotted line: 
Contribution from the discrete bound-state level, /"^(2;,Pt) 
Eq. (|4.2[) , shown with opposite sign. Solid line: Total distri- 
bution according to Eq. (|5.66[) . [Self-consistent soliton profile 
Eq. (53) with M = 0.35 GeV, M N = 3.26 M.] 



above p\ = 5M 2 and decreases rapidly at larger pt- 
This justifies our earlier use of the gradient expansion to 
study the large-pT behavior. At low pt there are very 
significant cancellations between the gradient expansion 
and the discrete level contribution, causing the sum to 
be 3-4 times smaller than the individual terms. Since 
the gradient expansion provides only a rough approxima- 
tion to the Dirac continuum at low px (cf. Fig. [FJ which 
compares different variants), and is subject to some un- 
certainty resulting from the UV cutoff, we cannot use 
Eq. ([51)61) t0 predict the total with any meaningful 
relative accuracy at low px, but can conclude only that 
it is substantially smaller than the gradient expansion 
result alone. More quantitatively, if we require that the 
level contribution be < 50% of the gradient expansion we 
conclude that we can safely use Eq. (15.661) for a numerical 
estimate at p\ > 2 M 2 . At lower values of pt one should 
use methods based on exact summation over levels to cal- 
culate the sea quark transverse momentum distribution 
[40l | . Note, however, that the contribution of this region 
to the pT-integral determining the total sea quark den- 
sity is rather small (see below), so that the question of 
the exact behavior of the sea quark p? distribution in 
this region is somewhat academic. 
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FIG. 14. Transverse momentum distributions of flavor-singlet 
unpolarized valence and sea quarks at x = 0.1. Panel (a) 
shows fu+<i-u.-d and ju+d ag f unct j ons f p 2, on a logarithmic 

scale; panel (b) shows the radial distribution 2-KpTfi +d ~ u ~ d 
and 2npTfi +d on a linear scale, such that the area un- 
der the curves corresponds to the integral over pr- Dashed 
lines: Valence quark distribution (see Fig.|6j). Solid 
lines: Sea quark distribution (PV regularization). [Self- 

consistent soliton profile Eq. (|A4() with M = 0.35 GeV, M N = 
3.26 M.l 



I. Sea vs. valence quark distribution 



Using the numerical approximation of Sec. IV Hi we now 
want to compare our results for the sea quark transverse 
momentum distribution with those of the valence quarks 
calculated in Sec. IIVI Figure [TJ] summarizes the numer- 
ical results for the valence distribution f^ +d ^ u ^ d (x,pr) 
and the sea quark distribution /" +c! (x,pt) at a represen- 
tative value of x = 0.1. Panel (a) shows the distributions 



themselves on a logarithmic scale; panel (b) the radial 
distributions on a linear scale, such that the area un- 
der the curves corresponds directly to their integral over 
Pt- Similar results are obtained at other values of x: 
the shape of the individual px distribution changes little 
with x (cf. Fig. U for the valence distribution); only their 
normalization changes in proportion to the total valence 
and sea quark density. 

The numerical estimates clearly show very different 
shapes of the valence and sea quark transverse momen- 
tum distributions, especially at large values of pr, as 
first observed in the calculation of Rcf. [4(i| . Based on 
our theoretical analysis we can now explain this strik- 
ing behavior as the effect of dynamical chiral symmetry 
breaking in the QCD vacuum on the intrinsic transverse 
momentum distribution of the sea quarks. Even with the 
strong modification of the would-be \jp\ tail by the UV 
cutoff, the sea quark transverse momentum distribution 
in the chiral quark-soliton model is qualitatively differ- 
ent from that of the valence quarks. While the precise 
numerical values depend on the model implementation 
(see e.g. Fig. ITT)) , the fact as such is rooted in the basic 
structure of the effective dynamics chiral and should be 
model-independent. 

When interpreting the results of Figure [14] one should 
keep in mind that the accuracy of the approximation 
Eq. (|5.66p used in our numerical estimate of fi +d {x,pr) 
is not sufficient to predict the values at p\ < 2 M 2 
with meaningful relative accuracy (cf. the discussion in 
Sec. IV Hp . In this sense the plot of the radial distribu- 
tion, in which the low-p^ region is suppressed, conveys a 
more realistic picture. This uncertainty, however, in no 
way influences our conclusions regarding the qualitatively 
different behavior of valence and sea quark distributions 
at large pr- 

The qualitative difference between the pt distribution 
of valence and sea quarks is the most important practical 
result of our study. Its numerous implications for deep- 
inelastic processes are explored in Sec. I Villi 



J. Polarized sea quark distribution 

To complete our study of the sea quark transverse 
momentum distribution we want to investigate also the 
flavor-nonsinglet polarized sea quark distribution. The 
gradient expansion of this distribution can be carried out 
in complete analogy to the flavor-singlet unpolarized case 
starting from Eq. (pOSf . cf. Sees. IV Al and lVBI we do not 
present the intermediate steps here. The result can again 
be represented as a convolution integral over the momen- 
tum of the classical chiral field, analogous to Eq. (|5.16l) . 



dy 



d 2 k T g c i(y, k T ) 



9qg(x,y;PT,k T ). 



(5.67) 
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The relevant momentum distribution of the classical field 
is now 



9c\{y,k T ) 



^_^tr fl [r 3 C/ cl (fe)C/ cl (fc)t] (5.68) 
[k = (k T ,yM N )}. 



The momentum distribution resulting from the quark 
loop integral turns out to be the same as in the flavor- 
singlet unpolarized case 



9qq{x,y;p T ,k T ) = f q q(x,y;p T ,k T ) 



(5.69) 



(for the jPr~integrated distributions this was already 
noted in Ref. [32j). This remarkable fact can be un- 
derstood as an instance of "restoration of chiral sym- 
metry." When expanding in gradients of the chiral fields 
the numerator of the quark loop integral becomes inde- 
pendent of the dynamical quark mass M, and the same 
coefficient is obtained for the axial vector-type operator 
in the polarized distribution (Dirac matrix 7 + 75) as for 
the vector-type operator in the unpolarized distribution 
(7 + ). A more microscopic interpretation of this result in 
terms of pair correlations in the nucleon wave function 
will be discussed below. Equation (|5.69[) has important 
consequences for the behavior of the flavor-nonsinglct 
polarized pt distribution at large transverse momenta. 

At transverse momenta p\ ^> M 2 the convolution in- 
tegral Eq. (|5.68[) can again be simplified by using the 
collinear approximation, Eqs. (|5.19p and (|5.20[) , and be- 
comes 

9i~J 3 a{x,Pt)= ( — 9d(y) g q q(x/y,p T ). (5.70) 

j x v 

Here g c \{y) is the corresponding fc<r-integrated light-cone 
momentum distribution of the classical chiral field, cf. 
Eq. ([Oil, 



9c\{y) 



d 2 k T g c \(y, k T ) 



T tr fl [r 3 [/ cl (fc)C/ci(fc) t ] (5.71) 



3(2tt) J (2tt) 
(fc 3 = yM N ), 

which can be evaluated using Eq. (|C6|) in Appendix [C] 
The numerical distribution is shown in Fig. [9l One sees 
that at y — >• the "polarized" distribution is suppressed 
relative to the "unpolarized" one, as expected, and that 
at large values of y ~ 1 the distributions approximately 
satisfy 3g c \(y) m f c \(y), corresponding to saturation of 
the large-iV c inequality for the transverse momentum 
densities (see below). The momentum distribution of the 
quark-antiquark pair appearing in Eq. (|5.70|) is the same 
as in the unpolarized case, cf. Eq. (|5.69[) , 



gqq{z,p T ) = fqq{z,p T ), 



(5.72) 



where z = x/y and f q q(z,pT) is defined in Eq. (|5.23|) and 
explicitly given by Eq. (|5.36l) , up to modifications by the 
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FIG. 15. Transverse momentum distributions of flavor- 
nonsinglet polarized valence and sea quarks at x = 0.1. 
Panel (a) shows g™- d - u + d anc i g^ d as functions of p% on 
a logarithmic scale; panel (b) shows the radial distribution 



and 27rpT3i on a linear scale, such that 



the area under the curves corresponds to the integral over 
Pt- Dashed lines: Valence quark distribution g^- d ~ u + d ( see 
Fig. [6]). Solid lines: Sea quark distribution g^ d (PV reg- 
ularization). [Self-consistent soliton profile Eq. (|A4|) with 
M = 0.35 GeV, M N = 3.26 M.] 



UV cutoff. As a result, the flavor-nonsinglct polarized 
distribution exhibits a would-be power-like tail at large 
transverse momenta similar to the the unpolarized dis- 
tribution, 



gt d (x,PT) 



ci- d (x) 



Pt 



[p 2 T > M 2 ). (5.73) 



Moreover, the coefficient of this tail is related to the 
UV divergence of the p^-intcgrated distribution in the 
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same way as in the unpolarized case, and is given by [cf. 
Eq. (533) ] 

Thus, our earlier discussion of the UV regularization and 
its effect on the transverse momentum distributions can 
be carried over directly to the polarized case. 

For a numerical estimate of the flavor-nonsinglet po- 
larized sea quark distribution at all values of px (includ- 
ing px ~ M) we use the analogue of the "interpolation 
formula" Eq. (|5.66j) . in which one adds the contribution 
from the discrete bound-state level, Eq. (|4.7p . to the gra- 
dient expansion approximating the Dirac continuum con- 
tribution, Eq. ([51)7]) . 

gt^PT) ~ 9t g L^PT) + 9%£(?,Pt)- (5-75) 

The resulting polarized sea quark distribution is shown 
in Fig. [15] and compared to the valence quark distribu- 
tion g^- d - a + d calculated in Sec. HVH One sees that, 
as in the flavor-singlet unpolarized case, the would-be 
power-like tail strongly influences the numerical behavior 
of the flavor-nonsinglet polarized sea quark distribution 
at p\ > few M 2 and causes it to be qualitatively different 
from that of the valence quarks. This again represents 
a direct effect of dynamical chiral symmetry breaking on 
the nucleon's partonic structure. Note that here the ef- 
fect occurs in a nonsinglet channel, in which the distribu- 
tions are likely to be much less affected by perturbative 
QCD evolution than in the singlet case, making this non- 
perturbative effect even more striking (cf. discussion in 
Sec. ED - 

A more microscopic explanation for the similarity of 
the flavor-singlet unpolarized and flavor-nonsinglet po- 
larized distributions at p\ ^> M 2 is provided in Sec. I VII 
where we show that the tails in the sea quark distribu- 
tions are due to correlated quark-antiquark pairs in the 
nucleon's light-cone wave function. The correlated pairs 
appear in scalar-isoscalar (E) and pscudoscalar-isovector 
(II) quantum numbers. The flavor-singlet unpolarized 
sea results from the overlap of like pairs (EE, nil), while 
the flavor-nonsinglet polarized one is due to the interfer- 
ence of different types of pairs (Ell, IIE) in the initial and 
final state. At p\ ^> M 2 the wave functions of the E and 
II pairs become the same due to the restoration of chiral 
symmetry, Eq. (|6. 2T[) , leading naturally to a relation be- 
tween the two sea quark distributions. We note that this 
derivation of the flavor-nonsinglet polarized sea quark 
distribution gives a precise meaning to the notion of "IIE 
interference," which was discussed in connection with the 
meson cloud model of flavor asymmetries in Rcf. [65]. 

To conclude our discussion of the polarized sea quark 
distribution, we would like to see how the general large- 
N c inequality for sea quark distributions, Eq. p.52[) . is 
realized in our model. At p\ ^> M 2 the distributions are 
given by the gradient expansion, in the simplified form 
of Eqs. (|5.21[) and (|5.70[) . and we want to verify that 

f? + *(x,pr) ± 3gt S (x,p T ) > (5.76) 



in this approximation. Because of Eq. (|5 . 72[) the differ- 
ences on the left-hand side of Eq. (|5.76[) can be written 
in the form 

/? +< Wt) ± -igt-\x,p T ) 
= f°- \Uv) ± 3 5cl (y)] f q ^x/y,p T ), (5.77) 

Jx v 

where fg q is explicitly positive; cf. the discussions in 
Sees . I V Dl and IVFl and the representation of this function 
as light-cone wave function overlap derived in Sec. IVI CI 
Eq. (|6.35p . The difference of the momentum distributions 
of the classical field in Eq. (|5.77[) . in turn, is given by 

F 2 M 2 N y f d 2 k T 

x tr fl [(l ±T 3 )U cl (k)U cl (ky\ (5.78) 
(fc 3 = yM N ), 

which is explicitly positive, cf. Eq. (|CJ7|) in Appendix [Ul 
Thus we see that the "restoration of chiral symmetry" ex- 
pressed in Eq. (|5.69[) naturally guarantees that the large- 
iV c inequalities for the sea quark distributions are satis- 
fied at p\ 3> M 2 in our scheme of approximations. 

The results of the interpolation formulas Eq. (|5.66[) and 
(|5.75[) for the distributions at p\ ~ AI 2 , taken literally, 
would violate the inequality Eq. (|5.76p at p\, < 2 M 2 . 
However, we noted in Sec. IV HI that in this region the 
unpolarized sea quark distribution cannot be estimated 
with any meaningful relative accuracy using this approx- 
imation. Wc therefore cannot conclusively study the in- 
equality at low px using this approximation. 



VI. SHORT-RANGE CORRELATIONS OF 
PARTONS 

A. Nucleon wave function at large momenta 

So far we studied the properties of valence and sea 
quarks in the nucleon by investigating their one-body 
momentum densities, Eqs. (|3.2[) and p.3[) . A more mi- 
croscopic understanding of our results can be obtained by 
considering the nucleon's partonic (or light-front) wave 
function in the chiral quark-soliton model. Specifically, 
we want to show that in this model sea quarks partly ex- 
ist in correlated pairs of a size of the order of the cutoff 
scale A -1 <C R, and that the "tail" in their transverse 
momentum density can directly be attributed to these 
configurations. That is, the two-scale picture of the ef- 
fective dynamics resulting from chiral symmetry break- 
ing implies the existence of short-range quark-antiquark 
correlations in the partonic wave function. This observa- 
tion has far-reaching implications for our understanding 
of the partonic structure of the nucleon not only in the 
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chiral quark-soliton model but in QCD in general. It sug- 
gests an interesting analogy with short-range nucleon- 
nucleon correlations in nuclei, which give rise to high- 
momentum components of the nuclear spectral function 
governing single-particle knockout reactions and can be 
probed directly in multiparticle correlation experiments 
[44H46j | . A detailed treatment of parton short-range cor- 
relations due to dynamical chiral symmetry breaking and 
their implications will be the subject of a subsequent 
publication. Here we want to discuss only the aspects 
relevant to understanding the behavior of the intrinsic 
transverse momentum distributions. 

The light-front wave function of the nucleon in the chi- 
ral quark-soliton model was derived and discussed in a 
general context in Refs. [U [42|. In this approach the 
many-body wave function of the fast-moving nucleon 
is constructed by applying the creation operators of N c 
valence quarks and a coherent superposition of quark- 
antiquark pairs to the vacuum state of the effective chi- 
ral theory, i.e., the Dirac vacuum with baryon number 
zero. The construction can be carried out explicitly in 
the sense of a Fock state expansion, usin g th e bound- 
state level occupied by N c quarks as a seed [43[ . Here we 
are interested in the sea quark component of the nucleon 
wave function at large transverse momenta, 



p 2 T > M 2 . 



.1) 



This component of the wave function can be calculated 
directly in a simple approximation, assuming the exis- 
tence of a stable mean field. It will be seen that this 
approximation is equivalent to the gradient expansion 
for the one-body momentum densities of sea quarks de- 
scribed in Sec. El 

The appearance of sea quarks in the light-cone wave 
function of the nucleon in the chiral quark-soliton model 
can be viewed as the creation of quark-antiquark pairs 
by the chiral field produced by the other constituents. 
In general this is a complicated process, which affects 
the state of motion of the source particles by chang- 
ing their longitudinal and transverse momenta. Also, 
the produced pair can re-interact with the chiral field 
and experience distortion of its wave function. The limit 
of large transverse momenta, Eq. (|6.1|) . permits several 
important simplifications, which make it possible to de- 
scribe this process in practice and calculate the sea quark 
component of the wave function in a controlled approxi- 
mation. 

First, the invariant mass difference in the pair creation 
process is dominated by the invariant mass of the pair, 
while the contribution from the change of the state of 
motion of the source system can be neglected. Consider 
the creation of a pair with overall plus momentum frac- 
tion y and transverse momentum kx by a "source" with 
mass ~ Mm (the precise coefficient does not matter) and 
initial transverse momentum zero. Let z and z be the 
relative plus momentum fractions of the quark and anti- 
quark in the pair, and pt the transverse momentum of 
the quark (cf. Sec. IV Bl) . The change of invariant mass 



of the total system (source and quark-antiquark pair) in 
the process is 



Astot 



M 2 +p 2 T M 2 + (k T -p T ) 2 



yz 

M 2 N + k 2 , 



M 



V 



N 



M 2 + p\ M 2 + (kr ~ Pt? 



y 2 M 2 N + yk 2 T 

y 



(6.2) 



The longitudinal and transverse momentum transfer by 
the source is of the order 



|fcr| 

yM N 



M, 
M. 



(6.3) 
(6.4) 



In the limit pf, ^> M 2 we can neglect the last term in 
the parenthesis in Eq. (|6.2p and drop kx and M in the 
other terms. The total invariant mass difference is thus 
determined by the intrinsic invariant mass of the quark- 
antiquark pair, s, cf. Eq. (|5.33[) . and given by 



As to 



Pt 

yzz 



(6.5) 



Second, in the region p\ ^> M 2 the interaction of the 
quark-antiquark pair with the source is effectively given 
by the leading-order Born approximation in the nontriv- 
ial part of chiral field, 



M[U l5 (x) - 1], 



(6.6) 



describing the deviation from the vacuum. Higher-order 
interactions come with higher powers of the constituent 
quark mass (as well as gradients of the chiral fields) and 
are suppressed by inverse powers of px- This means that 
at p\ M 2 the distortion of the internal wave function 
of the quark-antiquark pair can be neglected, and that 
it can be regarded as being in a plane-wave state after 
its creation. 

Altogether, we sec that a very simple picture of the sea 
quark component of the nucleon's partonic wave function 
emerges at p\ 3> M . The fast-moving nucleon creates 
color-singlet quark-antiquark pairs through the interac- 
tion Hamiltonian 

H int (t) = [ d 3 x$(t,x)M[U^{t,x)-l] ret iP(t,x), 



(6.7) 

where [U* 15 — l] re t is the nontrivial chiral field produced 
by the other constituents, in the sense of a retarded po- 
tential. The interaction is to be treated in first order, 
and the pair is produced in a plane-wave state. In this 
approximation the rest of the nucleon producing the pair 
acts as a classical source (henceforth called "classical nu- 
cleon" ) ; it transfers longitudinal and transverse momen- 
tum as well as spin/isospin quantum numbers to the pair 
but otherwise remains inert. 
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FIG. 16. (Color online) Approximation to the nucleon's light- 
cone wave function at pip 3> M 2 . The nucleon acts as a clas- 
sical source, which produces a color-singlet quark-ant iquark 
pair through the retarded potential Eq. (|6.7p . The longitudi- 
nal momentum fractions and transverse momenta in the final 
state are indicated on the right (y = 1 — y etc.). The out- 
going classical nucleon can be in a rotational state with N 
or A quantum numbers. The quark-antiquark pair can have 
scalar-isoscalar (£) or pseudoscalar-isovector (II) quantum 
numbers. 



One may ask whether the approximations proposed 
here respect chiral invariance. As it stands, first-order 
Born approximation with the Hamiltonian Eq. (|6.7[) docs 
not result in chirally invariant amplitudes. However, 
below wc shall see that chiral invariance is effectively 
restored in the momentum density of quarks and anti- 
quarks at high px, as it should be. 

With the above approximations it is straightforward 
to calculate the high-p^ component of the nucleon's par- 
tonic wave function (sec Fig. I16|) . Using ordinary time- 
ordered perturbation theory, we calculate the amplitude 
for a fast-moving nucleon at t = — oo to produce a quark- 
antiquark pair at time t = via the interaction Eq. (|6.7p , 
assuming that the interaction is adiabatically switched on 
starting from t = — oo. We obtain 

(q(pi,(Ti,fx,Ci) q(p 2 , cr 2 , f 2 , c 2 ), R(Pr, Ir) | 
o -| 

N(P N ,I N ) 



(-») / dtH int (t) 

J —OO 

{2nY S i3 \ Pl +p 2 + P R -P N )d, 

u(pi,(Tl,fl) T v(p 2 ,(T2,f2) 



Ex + E2 + Ef, 



E 



N 



Here the momentum and energy of the nucleon in the 
initial state arc 



Pn 

En 



(0, P), 



P 2 



M 2 



(6.9) 
(6.10) 



and the momenta and energies of the quark and antiquark 



Pi = (pr, yzP), 
pi = (k T -pt, yzP), 



E 



1,2 



P\,2\ 



M 2 , 



(6.11) 
(6.12) 

(6.13) 



such that the total momentum of the pair is 

k = P1+P2 = (fcr, yP). (6.14) 

The momentum of the recoiling classical nucleon state 
R is fixed by the momentum-conserving delta function, 
and its energy is 



E T 



P-pi-p 2 \ 2 



M 2 
m N . 



3.15) 



The nucleon and quark momentum states are normalized 
according to the relativistic convention, 

(N(P')\N(P)) =2E n (2tt) 3 S^iP' -P), etc. (6.16) 



Eq. (|6.8p and all subsequent expressions are to be con- 
sidered in the limit P — > 00, where the matrix ele- 
ment represents the nucleon wave function in the infinite 
momentum frame, in the approximation specified above. 

In Eq. (|6.8[) , T denotes the transition matrix element of 
the retarded potential between classical nucleon states, 

r = (r(p r ,i r )\ [u^(o) - i] ret \n(p n ,i n )) 

= F- 1 (E + l7 5T a n a ), (6.17) 

where S and n a are the effective scalar and pscudoscalar 
field of the transition. The spin-isospin quantum num- 
bers of the states are summarily denoted by 7/v = 
(Sn,Tn,Tn3, SW3), where Sn = Tjv is the spin/isospin 
and T/V3 and Sn3) their projections, and similarly for 
Ir. Note that the recoiling classical nucleon can be in a 
different rotational state from the initial one; in the case 
of a transition induced by the quark-antiquark operator 
considered here it can have N or A quantum numbers 
(we explicitly evaluate the matrix element below). Fur- 
thermore, in Eq. (|6.8j) cri l2 and fi l2 are the spin and fla- 
vor quantum numbers of the quark and antiquark states, 
and u and v their plane-wave Dirac spinor and isospinor 
wave functions, normalized as uu = —vv = 2M. The 
color quantum numbers are denoted by ci^, and the pair 
is in a color-singlet state. 

To evaluate Eq. (|6.8j) further, we note that the energy 
denominator in the limit P — > 00 is proportional to the 
invariant mass difference between the initial and final 
state, Eq. 



Ei + E2 + Ef 



E 



N 



2P 



2yP' 



3.18) 



where we have used the approximation Eq. (|6.5[) in the 
last step. We see that we can write the right-hand side 
of Eq. $£E§ as 

u(pi,<Ji,fi) r ^(p2,g 2 ,/ 2 ) 

Ei + E2 + Er — En 
= %yP [S tps(z,p T ;(Tl,fi;(T2,f2) 

+ U a iPf l (z,p T ;aiJ 1 ;a 2 J 2 )]. (6.19) 

Here tps and tp^ are the infinite-momcntum-frame wave 
functions of a quark-antiquark pair with scalar-isoscalar 
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and pscudoscalar-isovector quantum numbers in the ef- 
fective chiral model, 



M u{p 1 ,a 1 )v{p 2 ,a 2 ) 
= ; o hf , 



F w s 



(6.20) 
(r a ) hh - (6-21) 



They are independent of the overall longitudinal momen- 
tum and depend only on the quark momentum fraction z 
and the transverse momentum p T (we neglect the trans- 
verse momentum of the pair as a whole, kx, relative to 
the intrinsic transverse momentum pr), as well as on the 
spin and flavor quantum numbers. 



B. Spin structure of pair correlations 

The explicit expressions for the plane-wave Dirac 
spinors of the quark and antiquark, with the spin pro- 
jections (Ti and a 2 defined relative to the (fixed) 3-axis, 
are 



u(pi,<7i) = 



v(p2,cr 2 ) 



/ \/ Ei + M w(ai) 
\ \Pi\ 



-VE 2 



-\JE 2 



_ \Pi\ 

M a 2 w(a 2 ) 



(6.22) 



3.23) 



where <j % (i = 1,2, 3) are the Pauli spin matrices (to avoid 
confusion with the spin quantum numbers we distinguish 
them by the hat) and w{ai^ 2 ) is a two-spinor with 



iu(cti, 2 = 1/2) 



-1/2) 



(6.24) 

Substituting these expressions and evaluating the bilinear 
forms in the limit P — > oo, we obtain the spin structure 
of the pair wave functions, Eqs. (|6.20|) and (|6.21[) . as 



ij'i 



M 

= 771 7== d /l/2 

b^yzz s 

x [(&p T )a- 2 + (2z - l)Ma 3 a% ia2 , (6.25) 



M 



(r a )f lh 

,"3~2 i ji/fi.2i 



F^y/ZZ S 

x £[-(<7p T )o- 3 o- 2 + Ma\ ±CT2 . 



3.26) 



In the terms proportional to pr the spin projections of 
the quark and antiquark are parallel (the matrix is di- 
agonal); they correspond to configurations with orbital 
angular momentum L = 1. The presence of these con- 
figurations is a direct consequence of the chirally-odd 
structure of the coupling of the quark-ant iquark pair to 
the chiral field (1,75). Note that these terms have the 
same coefficient in both pairs, up to a phase factor. In 



the terms proportional to M in Eqs. (|6.25j) and (|6.26j) the 
spin projections of the quark and antiquark are antipar- 
allel (the matrix is off-diagonal); in the case of 4>n this 
term would correspond to the spin-flavor wave function 
of the pion in the nonrelativistic quark model. 



C. Restoration of chiral symmetry 

An interesting simplification happens with the pair 
wave functions of Eqs. (|6.25[) and (|6.26p in the region 
\Pt\ ^ M, where we want to use them to in our approx- 
imation to the nucleon wave function. In this region the 
first term in the spin wave function of Eqs. (|6.25|) and 
(16.261) dominates, and one finds 



4MX ,2 
Flzzs* ~ IVpair| 



3.27) 



(here a = 1,2,3 is fixed; no summation over a). Thus, 
chiral symmetry is effectively "restored" in the quark- 
antiquark pair wave function at high momenta. This 
finding is important in ensuring the chiral invariance of 
the one-body density of sea quarks, as shown below. It 
also explains the close connection between the flavor- 
singlet unpolarized and the flavor-nonsinglet polarized 
sea quark distributions observed in Sec. IV J I 

Equation (|6.27|) also shows that the high-p^ behavior 
of the sea quark density is entirely governed by the L = 1 
component of the pair wave function induced by chiral 
symmetry breaking. This fact has implications for our 
general understanding of chiral symmetry breaking in a 
light-front wave function description of the nucleon. 



D. Momentum distribution from pairs 

Let us now calculate the one-body momentum den- 
sity of sea quarks from the first-order "wave function" of 
Eq. fl6J§. Generalizing Eqs. (|3~2|) and (f373|) . the flavor- 
singlet unpolarized antiquark density is obtained as 



P 



(N(P',I 



N 



(2tt) 3 

: V Yl b* f „(p)b fa (p)\N(P t I N )) 



E 

f=u,a 



= 2E N (2n) 3 S^(P' - P) f? +i (x,p T ) (6.28) 

[p = (p T ,xP)] . 

The delta function on the right-hand side appears be- 
cause, in difference from Eqs. (|3.2[) and Q3.3p . the center- 
of-mass motion of the nucleon is now quantized and the 
states normalized according to Eq. (|6.16|) . Substituting 
for the nucleon states the first-order quark-antiquark 
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component of Eq. (|6.8[) . and resolving the constraints re- 
sulting from momentum conservation, we get 



f! +d (x,PT) 



P 



d 3 k (2yP) 2 N C A 
(2tt) 3 J (2tt) 3 2Ei2E 2 2Er2E 



N 



dy Id kx 



(2ir) & Ay 2 zz 



3.29) 



(6.30) 



where 



^EE E 



Ir fi,fa fijOU 



(6.31) 



In Eq. (|6.29p we have chosen the pair momentum fc as 
integration variable, cf. Eq. (|6 . 14[) . so that y = k 3 /P 
and z = x/y under the integral; in Eq. (|6.30[) we have 
changed the integration variable from k 3 to y and re- 
placed the energies by their values in the P — > oo limit, 
with En = (1 — y)P « P in our approximation. The 
factor N c results from the summation over the quark- 
ant iquark colors. The sum over spins and flavors in 
Eq. (|6.3ip is readily evaluated. Because of the differ- 
ent symmetry of the II and S pair wave functions the 
cross terms vanish, and we have 



.4 



EEE 



: E( S * E 

Ir 



n a *n Q )|Virl 2 (pt^m 2 ). (6.32) 



In the last step we have taken into account that in the re- 
gion p|, 3> M 2 the pair wave functions become alike due 
to the "restoration of chiral symmetry," cf. Eq. (|6.27[) . 
As a result, the transition II and £ fields effectively ap- 
pear in a chirally invariant combination, rendering our 
result chirally invariant. 

Combining everything, we can represent Eq. (|6.30l) in 
the form of a convolution integral of momentum densi- 
ties, as obtained earlier from the gradient expansion, cf. 
Eq. (IQT1) . 

ff +S {x,p T ) = f cl (y) f^(x/y,p T ), (6.33) 

j x v 



where now 



/ci(y) 



fqq{z,PT) 



2(2tt) 3 



/ d 2 k T ^(E*E + n a *n a ), 

Ir 



pair | 



(27r) 3 2zz 



(z = x/y). 



(6.34) 
(6.35) 



These formulas establish the desired connection with the 
gradient expansion of the one-body densities in Sec. [V] 
The light-cone momentum density of the chiral field 
Eq. (|5.22[) appears as the square of the effective scalar and 



pscudoscalar fields associated with the transition between 
the classical nucleon states in which the quark-antiquark 
pair was produced. The momentum distribution of the 
quark-antiquark pair, Eq. (|5.23p . in turn, is given by 
the squared light-cone wave function of the pair, accord- 
ing to the standard overlap formula (the particular factor 
1/[(2tt) 3 2zz] is a consequence of our definition of the wave 
function) . 

The UV cutoff of the effective dynamics can be in- 
troduced into our treatment of pair correlated pairs in 
the nucleon's light-cone wave function in the same man- 
ner as discussed in Sec. IV Fl The invariant mass cutoff 
Eq. (|5.56p can naturally implemented at the wave func- 
tion level by multiplying the pair wave functions with the 
cutoff function, 



(6.36) 
(6.37) 



In fact, the physical significance of this regularization 
scheme emerges here at the wave function level. Note 
that the restoration of chiral symmetry at 3> M 2 is 
preserved if the S- and II— type pairs are suppressed at 
large s in the same manner. 

It remains to show that the light-cone momentum den- 
sity defined by Eq. (|6.34[) coincides with the expression 
Eq. (|5 . 22[) obtained from the gradient expansion of the 
one-body density. Here we can use the fact that the 
transition matrix element defined by Eq. (|6.17p is invari- 
ant under longitudinal boosts and behaves like a light- 
cone wave function. The classical fields are (pseudo-) 
scalars, and the only effect of a longitudinal boost is on 
their space-time dependence, while the transition ma- 
trix element depends only on the four-momentum trans- 
fer between the initial and final classical nucleon state, 
k = P/v — Pr- More explicitly, the boost invariance of 
Eq. (|6.34p can be demonstrated by considering the case 
of small fields that are generated perturbatively from a 
pointlike classical source via a Lorentz-invariant retarded 
Green function. 

We can therefore evaluate Eq. (|6.34p in the rest frame, 
where the initial and final baryon 3-momenta are 



Pn 
Pr 



0, 



-yM N ), 



(6.38) 
(6.39) 



and the 3-momcntum transfer to the nucleon is 



k = P N P R = (k T ,yM N ). (6.40) 

In the sense of the 1 /N c expansion its components are of 
the order 



yM N , \kr\ = 0(N°), 



3.41) 



and the energy difference between the initial and final 
classical nucleon states is 



En - M 



N 



k\ 2 
k\ 2 



2M N 



M 2 N - M N 



OiN- 1 ). 



3.42) 
3.43) 
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In this frame the classical nucleon effectively behaves as 
a nonrelativistic system. In particular, as a consequence 
of Eq. (|6.43|) , the "plus" momentum transfer between the 
states is in leading order of 1/N C expansion completely 
given by the 3-momentum transfer, and we obtain 

k+ = P+ - R+ = M N - E R - R 3 
= yM N + OiN- 1 ) 

= yP+ + OiN' 1 ), (6.44) 

as it should be, cf. the discussion in Sec. IIII CI In this 
frame the space-time dependent chiral field is given by 

lP'(t,x) nt = U%(x-X), (6.45) 

where X is the center-of-mass coordinate of the classi- 
cal nucleon, and the matrix element between momentum 
eigenstates of the classical nucleon is calculated as [3(| 

(p r \[u^(o)-iu\p n ) 

= 2M N J d 3 X e KPR-PN)X pg( x _ x) - 1] 



2M N Wj{k) k=(yMN M 
'1 + 75. 



2M 



N 



-U cl (k) 



1 



75 



(6.46) 



2 " w 2 

cf. Eq. (|5.6|) : the factor 2Mm results from the relativistic 
normalization of states, Eq. (|6.16[) . Finally, the matrix 
element between rotational states of the classical nucleon 
is obtained by subjecting the classical field to an (iso-) 
spin rotation 

U cl ->■ RU cl R) (6.47) 

and calculating the transition matrix elements between 
the initial and recoiling rotational states, using 



In) = 



dR(f> lR (R)* ... 4>i N (R), (6.48) 



where the integration is over the group measure the rota- 
tional wave functions are given in terms of the the Wigncr 
D-functions [30| • In the case at hand we actually do not 
need to calculate the integral Eq. (|6.48[) for a given Ir, as 
the sum over rotational quantum numbers in Eq. (|6.34[) 
produces a delta function that rigidly couples the rota- 
tions in the matrix clement and its complex conjugate, 



= S(R-R') 



(6.49) 
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The rotational average thus reduces to a single integral as 
in the expectation value of the one-body density. Com- 
bining everything, we finally obtain 



Ir 

F= 



E*E + n a *n a ) 
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tr [r f r] 



2M 2 N F 2 tr fi [U cl (k)U cl (kY} 



(6.50) 



Inserting this result into Eq. (|6.34[) we obtain precisely 
Eq. (|5.22p . This completes the proof that the first-order 
nucleon wave function Eq. (|6.8j) reproduces the gradient 
expansion result for the one-body momentum density of 
sea quarks discussed in Sec. |V] 

The calculation of the sea quark density from cor- 
related pairs in the nucleon's light-cone wave function 
presented here can easily be extended to the flavor- 
nonsinglct polarized distribution. One finds that this 
distribution originates from the interference of E- and 
H-type pairs in the wave function of the initial and final 
state; this structure was already discussed in connection 
with the meson cloud model in Ref. J65|. The restora- 
tion of chiral symmetry in the light-cone wave function 
of the pairs [cf. Eq. (|6.27[) for the probabilities] natu- 
rally leads to the result that the £11 and III] wave func- 
tion overlap at p\ ^ M 2 is of the same form as the 
EE and nil ones producing the flavor-singlet unpolar- 
ized density. In the conventions of Sec. IV Jl this implies 
that g q q{z,pT) = fqq(z,p T ), Eq. (|5.72p . which was ob- 
tained independently from the gradient expansion of the 
respective densities. Furthermore, the pertinent transi- 
tion fields E and II a , now defined as matrix elements in 
which the initial/final nucleon state has a definite spin 
projection 1/2 and isospin projection 1/2, combine to 
produce the polarized momentum distribution Eq. (|5.68p . 
Thus, the existence of E and II pair correlations and the 
restoration of chiral symmetry at p\ 3> M 2 naturally ex- 
plain the close connection between the flavor-singlet un- 
polarizcd and the flavor-nonsinglet polarized sea in the 
chiral quark-soliton model. 



VII. SUMMARY AND DISCUSSION 
A. Summary of model results 

In this article we have studied the effect of dynamical 
chiral symmetry breaking on the intrinsic transverse mo- 
mentum distributions of partons, using the chiral quark- 
soliton model as an approximate description of the effec- 
tive dynamics below the chiral symmetry-breaking scale. 
To conclude our investigation, we would like to summa- 
rize the findings from our model calculations, revisit the 
question of the embedding in QCD, and outline possible 
applications of our results to deep-inelastic processes. 

The results of our study of transverse momentum dis- 
tributions in the chiral quark-soliton model can be sum- 
marized as follows: 

• The constituent quark picture of the effective chi- 
ral dynamics implies a natural definition of the in- 
trinsic transverse momentum distributions of par- 
tons, extending the established description of pt~ 
integrated parton densities in this approach. The 
model describes the pT~distributions of constituent 
quarks and antiquarks, which are to be matched 
with QCD quarks, antiquarks and gluons at the 
chiral symmetry-breaking scale. 
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The distribution of valence quarks (quarks minus 
antiquarks) is dominated by momenta of the or- 
der of the inverse nucleon size, pr ~ M ~ R^ 1 
and shows an approximate Gaussian shape. The 
systematic differences between the unpolarized and 
polarized valence quark pr distributions attest to 
the relativistic character of the nucleon bound 
state. 

The distribution of sea quarks is qualitatively dif- 
ferent and shows a would-be power-like tail ~ 1 jp T 
extending up to the UV cutoff. Its coefficient is de- 
termined by low-energy chiral dynamics and quasi 
model-independent. Such behavior is found in 
both the flavor-singlet unpolarized and the flavor- 
nonsinglct polarized distribution. 

The UV cutoff of the model influences the shape 
of the sea quark transverse momentum distribu- 
tion at momenta p T ~ fewM 2 . Imposing rather 
general physical conditions on the regularization 
scheme (analyticity, charge conservation, large- 
distance behavior) stable numerical results for the 
sea quark distributions are obtained up to p T ~ 
10 M 2 . The sea quark px distribution thus ob- 
tained is qualitatively different from that of the va- 
lence quarks and extends up to much larger values 
of pt- 

The power-like tail in the sea quark distribu- 
tion can be explained microscopically as the re- 
sult of short-range correlations in the nucleon's 
light-cone wave function. The relevant configu- 
rations are quark-antiquark pairs with transverse 
momenta p T 3> M 2 created by the classical chi- 
ral field. These pairs can have scalar-isoscalar 
(£) and pseudoscalar-isovector (IT) quantum num- 
bers, whose wave functions become identical at 
p T ~> M 2 , corresponding to an effective "restora- 
tion of chiral symmetry" at the cutoff scale. 

The coordinate space correlation function obtained 
as the Fourier-Bessel transform of the sea quark pt 
distribution is unambiguously defined in the model. 
At large transverse distances it decays exponen- 
tially with a characteristic range determined by the 
constituent quark mass and the spatial size of the 
mean field. At zero transverse distances it coin- 
cides with the pT-integrated parton density. It is 
thus tightly constrained and largely independent of 
the regularization scheme. 

The flavor-singlet unpolarized and the flavor nons- 
inglet polarized sea quark pt distributions at p T ^> 
M 2 satisfy the large- N c inequalities for transverse 
momentum distributions, as a result of the in- 
variance properties of the classical chiral field and 
the "restoration of chiral symmetry" in quark- 
antiquark pairs with p\ 3> M 2 . 



When interpreting the results of our study or using 
them for phenomenology, one should keep in mind that 
the specific numerical results for the valence and sea 
quark transverse momentum distributions (see Figs. [H] 
and I15[) are model-dependent and should not be taken 
too literally The numerical values of the sea quark trans- 
verse momentum distributions are strongly affected by 
the UV cutoff, and while it is very encouraging that the 
physically motivated regularization schemes described 
here produce stable distributions, we cannot exclude that 
further refinement of the UV regularization would change 
the numerical values. Furthermore, as explained in Sec.|I] 
and discussed further in the following subsections, the 
distributions refer to effective degrees of freedom defined 
by the model, and we presently cannot express them ob- 
jectively as matrix elements of a priori defined QCD op- 
erators in the nucleon. Our main result is the qualitative 
difference between the valence and sea quark transverse 
momentum distributions. We expect it to be model- 
independent because it is rooted in the basic structure 
of the effective dynamics resulting from chiral symmetry 
breaking and can at least qualitatively be explained in 
QCD proper (see Sec lVIICj) . Any applications of our re- 
sults should therefore focus on this qualitative difference 
rather than the specific numerical values of the distribu- 
tions (see Sec. lVIIIj) . 

The study of the pr distributions in the chiral quark- 
soliton model presented here could be extended in vari- 
ous directions. One obviously should calculate also the 
Pt distribution of the 1 /A^ c -suppressed combinations of 
the parton densities, namely the flavor-nonsinglet unpo- 
larized densities, and /f and the flavor-singlet 
polarized ones, <?" +d and These distributions are 
proportional to the angular velocity of the isorotational 
motion of the soliton, ft ~ 1/N C and the correspond- 
ing expressions can be derived in analogy to those of the 
leading distributions in Sec. IIIIB| using the techniques 
described in Refs. [32l [o6j . An interesting question is 
whether the 1 /N c -suppressed sea quark distributions also 
exhibit a power-like tail at p T ^> M 2 , and to what extent 
this behavior would differ from that of the l/Vc-leading 
distributions. It is known that the UV cutoff dependence 
of the flavor-nonsinglet unpolarized distribution in the 
model is qualitatively different from that of the singlet 
distribution at paramctrically small values x ~ M/(AN C ) 
[Hj], suggesting a qualitatively different behavior also of 
the pr distributions. 

Equally interesting is the pt dependence of the quark 
transversity distributions h\ and h\. Earlier studies 
fo7l HH revealed that the sea quark transversity distri- 
butions in the chiral quark-soliton model are finite in 
the limit of large UV cutoff and numerically small, which 
would indicate that there is no ~ tail in the pt dis- 
tribution. It would be especially interesting to explain 
the different behavior of the transversity distributions at 
the wave function level, using the method developed in 
Sec. ED 

Another direction is the study of short-range parton 
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FIG. 17. Unpolarized flavor-singlet parton momentum densities xf\ 



u-\-d — u—d 



(x) and xfi d (x). The sea quark momentum 



densities are plotted multiplied by a factor 2, with which they appear in the momentum sum rule, (a) Chiral quark-soliton 
model. Shown are the results obtained with the numerical approximations Eqs. (|4.3p and (|5.66|l [self-consistent soliton profile 
Eq. (|A4[) . M — 0.35 GeV, Mm = 3.26 M; sea quark distribution with PV regularization]. The constituent quarks and antiquarks 
carry the entire light-cone momentum of the nucleon in the model, (b) QCD parton densities from the global fit of Ref. [34|] 
(GJV07 LO, scale fj^o = 0.3 GeV 2 ). Gluons carry ~ 30% of the nucleon's light -cone momentum. 



correlations at the wave function level. In Sec. IVII we 
calculate the nucleon wave function at p\ ^> M 2 in the 
simplest possible approximation, which violates chiral in- 
variance (although it is effectively restored at high px)- 
It would be worth developing a manifestly chirally in- 
variant expansion scheme, which takes into account the 
distortion of the pair wave function by the classical field. 
This could be done using the induced vector field repre- 
sentation of the effective chiral dynamics, in which the 
original matrix field C/ 75 is absorbed by a chiral rotation 
of the massive fermion fields. Furthermore, one should 
further explore the connection of our approximation with 
the Fock space expansion of the large-iV c baryon wave 
function of \4^k- 

B. Matching with QCD 

The approach to partonic structure presented here is 
based on the idea of an effective description of QCD 
below the chiral symmetry-breaking scale. The chiral 
quark-soliton model of the nucleon describes the trans- 
verse momentum distribution and correlations of con- 
stituent quarks and antiquarks — effective degrees of 
freedom, which arc to be matched with QCD quarks, 
antiquarks and gluons at the chiral symmetry-breaking 
scale. As explained in Sec. HI this matching cannot be 
performed entirely on the basis of intrinsic properties of 
the effective chiral dynamics but requires additional in- 
formation about its embedding in QCD, either in the 
form of a microscopic derivation of the effective model or 
of phcnomcnological input. A detailed treatment of this 



problem is beyond the scope of the present study. Still 
we would like to comment on several aspects and draw 
some conclusions based on the structure of the results 
obtained within the effective model. 

To discuss the matching quantitatively it is instructive 
to compare the ^"integrated parton densities in the chi- 
ral quark-soliton model, which were studied in several 
earlier works [32l HH, [56[, with empirical parametriza- 
tions of the parton densities obtained from global fits 
(sec Fig. fl7|) . Fig. ITTk shows the unpolarized flavor- 
singlet momentum densities of valence and sea quarks 
in the model, xfi +d ~ a - d {x) and xff +S {x). For consis- 
tency we show here the results obtained with the nu- 
merical approximations employed in the present study 
of pr distributions, Eqs. (|4.3[) and (|5.66[) : these approx- 
imations reproduce the exact numerical results for the 
PT-integrated valence and sea quark densities [56| with 
a relative accuracy far better than 10% and 20%, respec- 
tively, in the region x < 0.5. Note that the momentum 
sum rule is satisfied within the model in leading order of 
the 1/N C expansion [32[, which is another testimony to 
the consistency of the relativistic mean-field approxima- 
tion. The sum rule is preserved by the PV regularization 
used here; see Ref. [56| for a detailed discussion. The dis- 
tributions of constituent quarks and antiquarks obtained 
in the model thus satisfy 

J dxx[ff +d - a - d (x) + 2/f+ J (x)] model = 1. (7.1) 

In Fig. 117k the sea quark distribution is plotted including 
the factor 2 with which it appears in Eq. (|7.1[) . FigurcfTTb 
shows the empirical paramctrizations of the QCD parton 
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momentum densities obtained in a recent global fit of DIS 
and other data [34| . For clarity we show here the leading- 
order (LO) distributions, which arc independent of the 
renormalization scheme and can readily be compared 
with the model distributions. The scale of these distribu- 
tions is ^lo = O-^ GeV 2 , corresponding approximately to 
the chiral-symmetry breaking scale p~ 2 w 0.4 GeV 2 iden- 
tified in the nonperturbative QCD approaches discussed 
in Sec. HI In the QCD parton densities the momentum is 
shared between quarks, antiquarks and gluons, and the 
momentum sum rule reads 



dx x [f{ 



u-\-d—u—d/ 



(x) + 2/f+ d (x) + g(x)] QC B = 1. 



_ (7-2) 
According to the fit of Ref. [34[ ~ 30% of the nucleon's 
momentum at the low scale is carried by gluons. Com- 
paring the model and the QCD parton distributions in 
Fig. [TTk and b two features stand out. First, the total 
momentum carried by sea quarks is roughly the same; 
however, the constituent sea quark distributions are no- 
ticeably harder (stronger at larger x) than those of the 
QCD antiquarks. This shows that the matching of the 
model sea quark distribution with QCD partons is gener- 
ally nontrivial and may produce and evolution-like effect. 
[Note that the model distributions should not be consid- 
ered as genuine predictions in the parametrically small 
region x ~ M 2 /(A 2 7V C ), where they are sensitive to the 
details of the UV cutoff (see Sec. |Vj) .] Second, the mo- 
mentum carried by the valence quarks is smaller in the 
case of the QCD distributions. This suggests that most 
of the QCD gluons should be "generated" from the model 
valence quark distributions in the matching process. 

In the simplest approximation one can identify the con- 
stituent quarks and antiquarks of the effective model with 
QCD quarks and antiquarks at the chiral symmetry- 
breaking scale and set the gluon distribution in the model 
to zero. For the p-r-integrated parton densities this ap- 
proximation can be formally justified in the instanton 
model of the QCD vacuum, where it appears as the 
leading-order approximation in the expansion of the in- 
stanton packing fraction. Fig. [T7] and the above discus- 
sion show that this approximation is reasonable (at non- 
exceptional values of x) but of limited accuracy. The 
limitations of this approximation should be kept in mind 
when interpreting our results for the px distributions. 

An important lesson from the fit shown in Fig. 117b is 
that the empirical gluon density at the low scale shows 
substantial strength at large values of x > 0.2 and has 
a shape comparable to the valence quark distribution. 
Such behavior is difficult to explain as the result of a 
composite structure of individual constituent quarks and 
antiquarks, which, like DGLAP evolution, would pro- 
duce gluons primarily at values of x much smaller than 
that of the quarks. It suggests that this component of 
the gluon density originates rather from correlations be- 
tween constituent quarks, which can leverage the sum 
of the x-valucs of the two correlated quarks for the pro- 
duced gluons. This fact has not been appreciated in most 



attempts to explain the empirical gluon density made 
so far. Whether the relevant correlations are between 
quark-antiquark pairs in the found in the present 

model based on effective chiral dynamics, or diquark- 
like correlations between valence quarks, is an important 
question which deserves further study. 



C. Toward parton correlations in QCD 

The chiral quark-soliton model predicts short-range 
correlations between constituent quarks and antiquarks 
as a consequence of a basic property of the effective chi- 
ral dynamics: A 2 ^> Af 2 , or the parametric smallncss 
of the dynamical quark mass compared to the UV cut- 
off representing the chiral symmetry-breaking scale. An 
important question is whether and how such parton cor- 
relations could be rigorously defined in QCD. The chal- 
lenge lies in the fact that the relevant dynamical scale 
arises from chiral symmetry breaking in the QCD vac- 
uum, which is not readily associated with partonic struc- 
ture in a model-independent manner. Gribov's concept 
of the partonic wave function of a fast-moving hadron 
[69| maintains the connection between partons and vac- 
uum structure and should in principle be appropriate for 
discussing parton correlations as proposed here; however, 
the concept was developed based on scalar field theory, 
and the extension to QCD presents many technical chal- 
lenges (gauge invariance, UV divergences, renormaliza- 
tion) . A more rigorous formulation of parton short-range 
correlations in QCD may be possible with the new con- 
cept of multiparton distributions, which were introduced 
in the description of multiple hard scattering processes 
in high-energy pp collisions; their operator definition and 
renormalization prop erties have been studied in several 
recent works [7(| l7lf . 

Short of a rigorous formulation, we can still develop 
a qualitative picture of how parton short-range correla- 
tions emerge from dynamical chiral symmetry breaking 
in QCD. In the usual equal-time formulation of rclativis- 
tic dynamics the QCD vacuum is not empty, but popu- 
lated by localized nonperturbative gluon fields. These 
fields create quark-antiquark pairs with a characteris- 
tic size p ~ 0.3 fm [cf. Eq. (| 1 .4[) and following discus- 
sion], which form the chiral condensate (see Fig. [TSkV 
Quarks propagating through this medium interact with 
the vacuum fields and effectively acquire a dynamical 
mass, which determines much of hadronic structure. 
This phenomenon has been investigated extensively in 
slow-moving hadronic states (P < whose prop- 

erties can be studied using Euclidean (imaginary-time) 
correlation functions. A partonic description appears 
when considering hadrons which move with a momen- 
tum much larger than the scale of the vacuum fluctu- 
ations, P S> Following Gribov [6^|, in this limit 
one can separate the quanta carrying a finite fraction of 
the hadron momentum from those "left behind" in the 
vacuum, and the hadron becomes a closed system in the 
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FIG. 18. (Color online) Qualitative picture of the emergence 
of nonperturbative parton short-range correlations in QCD 
(details see text), (a) Slow-moving nucleon. (b) Fast-moving 
nucleon (parton picture). 



quantum-mechanical sense, amenable to a wave function 
description. When we imagine approaching the regime 
of large momentum P » gradually, it is clear that 
some of the quark-antiquark pairs in the vacuum will be 
"dragged along" and become the sea quarks in the nu- 
cleon's partonic wave function (see Fig.[T8b). These pairs 
of course inherit the typical size p with which they exist in 
the vacuum and thus induce nonperturbative short-range 
correlations in the nucleoli's partonic wave function. In 
this argument we implicitly assume that the nonpertur- 
bative wave function is defined "at a scale of the order 
p -2 ," and that configurations with transverse momenta 



Pt 



2 > 



will be built up by pcrturbativc QCD radia- 



tion. While leaving aside many important questions (UV 
divergences, renormalization) this simple picture qual- 
itatively explains how parton short-range correlations 
emerge from chiral symmetry breaking in QCD [72| . 

The picture described here implies that the small-size 
qq pairs in the nucleon wave function are generally ac- 
companied by strong gauge fields with transverse mo- 
menta of the order These gauge fields need not 
necessarily project onto physical gluon states in the limit 
P — > oo (e.g., they can correspond to unphysical polariza- 
tion states), but can exert forces on the quark and anti- 
quark corresponding to higher-twist effects. A transverse 
momentum-dependent hard scattering process involving 
sea quarks with px ~ p^ 1 thus generally takes place in 
the presence of a strong gluon field. One should therefore 
expect sizable corrections to the impulse approximation, 
in which one takes into account the px ~ p^ 1 of the 
initial quark/ antiquark but not the equally strong final- 
state interaction with the small-size gluon field. For this 



reason we cannot use our calculated px distributions di- 
rectly to make numerical predictions for the transverse 
momentum distributions in hard processes. Our conclu- 
sions regarding semi-inclusive DIS presented below as- 
sume only on the existence of the short-distance scale 
and do not rely on the impulse approximation. 

Modeling the effect of the nonperturbative gluon fields 
at the chiral symmetry-breaking scale in hard scattering 
processes will be essential for putting the proposed pic- 
ture of short-range correlations in QCD on a more quan- 
titative footing. The instanton vacuum model has proved 
to be a valuable tool for estimating the effect of such fields 
in inclusive DIS, where one can apply the local operator 
product expansion to identify the scaling (leading-twist) 
and power-suppressed (higher-twist) parts of the cross 
section [73Tj75|. Whether the instanton vacuum model 
could be adapted to estimate also final-state interactions 
in semi-inclusive processes is an interesting question for 
further study. 

Another interesting question concerns the color struc- 
ture of parton short-range correlations in QCD. In the 
effective model used in the present study the chiral 
symmetry-breaking color fields are integrated out, and 
the resulting dynamics generates short-range quark- 
antiquark correlations only in color-singlet states (£, II). 
The parameter N c in the model really plays the role of 
a degeneracy or weight factor and no longer refers to ac- 
tual color interactions. In a more microscopic approach 
such as the instanton vacuum, short-range correlations 
could in principle appear also between quark-antiquark 
pairs in color-octet states, if the color is compensated by 
physical gluon fields. Such pair correlations would have 
different properties as initial conditions for QCD radia- 
tion starting from the scale p~ 2 , and would influence the 
final state in semi-inclusive processes in a different way, 
compared to color-singlet correlations. Investigating the 
possibility of such color-octet correlations in the instan- 
ton vacuum model would be an interesting problem for 
further study. 

In sum, our arguments suggest that the short-distance 
scale associated with dynamical chiral symmetry break- 
ing in the QCD vacuum leaves a characteristic imprint 
on the nucleon's partonic structure. It implies a definite 
pattern of intrinsic transverse momentum distributions, 
in which the distribution of valence quarks is concen- 
trated at p T of the order of the inverse hadronic size 
R~ 2 , while that of sea quarks extends up to the chiral 
symmetry-breaking scale p~ 2 S> R~ 2 . At the level of the 
partonic wave function, the sea quarks show short-range 
pair correlations that reflect their creation by nonpertur- 
bative gluon fields of transverse size p. These conclusions 
are generic and rely only on the existence of the nonper- 
turbative short-distance scale in the QCD vacuum; they 
do not depend on the particular implementation of this 
scenario in the effective dynamical model employed here. 
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VIII. APPLICATIONS TO DEEP-INELASTIC 
PROCESSES 

Our findings about the intrinsic transverse momentum 
distributions of partons have many implications for DIS 
experiments with identified particles in the final state. 
A quantitative description of the transverse momentum 
spectra in particle production requires further assump- 
tions about perturbative QCD radiation and final-state 
interactions and detailed modeling of the fragmentation 
process, and is beyond the scope of the present article. In 
the following we only want to point out several obvious 
consequences and applications of our results that can be 
stated model-independently and merit further detailed 
study. 

A. Perturbative QCD radiation 

When applying our results to deep-inelastic processes 
with identified particles it is important to take into ac- 
count the effects of perturbative QCD radiation. The pic- 
ture of intrinsic transverse momentum distributions and 
parton short-range correlations in the nucleon wave func- 
tion described in Se dVIll applies at the chiral symmetry- 
breaking scale p~ 2 . Generally, in a hard process pertur- 
bative QCD radiation builds up configurations with in- 
variant masses (or virtualities) of the order of the hard 
scale, Q 2 . In inclusive DIS this radiation is described 
by DGLAP evolution and well understood. In processes 
where one measures the transverse momenta of particles 
in the final state the relevant radiation processes are gen- 
erally much more complex; for semi-inclusive DIS this 
problem was recently studied in Ref. (l4j . A key ques- 
tion is whether the nonperturbative short-distance scale 
implied by chiral symmetry breaking, p~ 1 « 0.6 GeV, 
acts as an infrared cutoff for perturbative radiation in 
such processes, or whether such radiation (possibly with 
Sudakov suppression) is still relevant at lower scales. Ex- 
perimental evidence on this point is ambiguous. Data 
on jet structure in e + e~ annihilation can be explained 
by allowing perturbative radiation to substantially lower 
scales. At the same time, there are hints from exclusive 
processes that low-virtuality radiation is suppressed, e.g. 
in the surprisingly similar magnitude and energy depen- 
dence of transverse and longitudinal cross sections in p° 
meson production [76|. What these observations imply 
for semi-inclusive particle production is an interesting 
question which merits further study. 

B. Semi— inclusive measurements 

Semi-inclusive DIS with single identified hadrons in 
the current fragmentation region, 7* + N — > h + X, is a 
standard tool for separating the different charge and fla- 
vor components of the nuclcon's parton densities, using 
the fact that the fragmentation process is sensitive to the 



charge and flavor of the struck quark. In these measure- 
ments one integrates over the transverse momentum Pt,k 
of the identified hadron h and aims to extract the cross 
section as a function of the fraction z% of the virtual pho- 
ton energy in the target rest frame carried away by the 
hadron. Generally, the observed transverse momentum 
PT,h is compounded from the intrinsic transverse momen- 
tum of the struck parton, QCD final-state interactions 
and perturbative radiation, and the transverse momen- 
tum incurred during the soft fragmentation process. Our 
findings about the intrinsic transverse momentum distri- 
butions imply that hadrons produced in scattering from 
antiquarks generally have a much broader Pt,h distribu- 
tion than those produced from quarks, if the transverse 
momenta incurred from final-state interactions and frag- 
mentation are comparable in both cases. One example 
is the production of K + (valence quark content us) and 
K~(us), where our picture predicts a broader Pt,h dis- 
tribution for the latter, assuming the production is dom- 
inated by scattering from u and u quarks in the proton. 
In experiments with incomplete coverage in Pt,k± this 
can result in a modification of the observed numbers of 
K + and K~ that is not related to the u and u num- 
ber densities in the target and must be corrected in the 
charge/flavor separation. 

Another instance where this effect plays a role are mea- 
surements where relations between fragmentation func- 
tions are used to isolate certain combinations of parton 
densities [z3, [l^- For example, the cross sections for 
semi-inclusive 7r + and 7r~ production, integrated over 
the pion transverse momentum Pt^-, is up to an overall 
kinematic factor given by 

+ e 2 [f^x)Di / ^(z 7T ) + f^x)Df / ^(z 7r )}, (8.1) 

where e u ^d are the quark charges and (z^), etc., the 

fragmentation functions describing the inclusive proba- 
bility of a quark/antiquark to produce a carrying 
fraction z^ of its longitudinal momentum, integrated over 
the soft transverse momenta. Taking the difference of ir + 
and ir~ cross sections, and using the relations between 
the fragmentation functions following from charge conju- 
gation invariance, 

D^ ± (z 7r ) = D^ T (z 7r ), 
D d ^ ± {z^)^D^\z^) 1 (8.2) 
and from isospin symmetry, 

D^ ± {z„)=D d x / ^{z„), 
D^ ± (z n )=Df / ^(z n ) ) (8.3) 

one obtains 

- « [el ~ e%m{x) - ff(x) - ff(x) + f?(x)} 
x [D™ /n+ (z„) - (z„)\. (8.4) 
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A measurement of this cross section difference thus pro- 
vides direct access to the valence quark densities in the 
target; the sea quark densities drop out because of the 
relations between the fragmentation functions. However, 
this reasoning requires modification if the effects of trans- 
verse momenta are taken into account and the experi- 
ment covers only part of the relevant Pr^-distribution. 
To illustrate the point, we may use the simple parton 
model with intrinsic transverse momenta in the parton 
density and the fragmentation function, where it is as- 
sumed that all transverse momentum integrals converge 
because of some intrinsic soft scale (our conclusion is 
model-independent and holds also in the presence of 
QCD radiation and final-state interactions). In this 
model the transverse momentum distribution of the iden- 
tified hadron is given by the convolution of the intrinsic 
transverse momentum distribution of quarks in the target 
and that incurred in the fragmentation process, 

u^{P T ,„) cx J d 2 PT J d 2 K T 6W(P T ,„ - z^ T - K T ) 

x [elf^x,p T )Dl /n± {z^K T ) + ...], (8.5) 

where the ellipsis denotes the corresponding other terms 

appearing in Eq. (|8. 1|) . Here (z^^Kt) etc. de- 

note the transverse momentum dependent fragmentation 
function, which satisfy 



jd 2 K T D^ + (z v ,K T ) = D^ + (z„), 



etc. (8.6) 



The i^T _ dependent fragmentation functions obviously 
obey the same charge conjugation and isospin symme- 
try relations as the ifr-integrated ones, Eqs. Q8.2[) and 
(|8.3[) . As a result, one obtains a formula analogous to 
Eq. (|8.5|) for the difference of cross sections measured at 
the same Pt,-k- 

[** + 

cx J d 2 PT Jd 2 K T 6^{P TtW - z^p T - K T ) 
x {el- e 2 d ) \K(x,p t ) - ff(x,p T ) 



ff(x, PT ) + tf{x iPT ) 



x [D u ^ (z„,K T ) 



D 



u/tt 



(z*,Kt)] 



.7) 



The cross section difference is proportional to the dif- 
ference of pr-dependent quark and antiquark distribu- 
tions. If one integrated over the pion transverse mo- 
mentum Pt.tt without restriction, the delta function in 
Eq. (|8.7j) would disappear, and the unrestricted integrals 
over px and Kt would reduce the transverse momentum- 
dependent parton distribution and fragmentation func- 
tions to the integrated functions of Eq. (|8.5p . However, 
one can no longer extract the px -integrated valence quark 
density from Eq. (|8.7j) from measurements with incom- 
plete Pt,tt coverage if quarks and antiquarks have dif- 



ferent intrinsic pt distributions, as implied by our ar- 
guments based on the QCD vacuum structure. To see 
this, let us write the pT _ dependent quark and antiquark 
distributions in the simple parton model in the form 

ff{x,p T ) = ft{x)F?{x,p T ), (8.8) 
f?(x,p T ) = f?(x)F?(x,p T ), (8.9) 

and similarly for the d quarks and antiquarks. Here the 
functions F" describe the normalized px profile of the 
quarks and antiquarks at a given x, 



r d 2 PT F?(x,pr) 



= 1, 



(a = u,u,d,d). (8.10) 



We can then write the differences of py-dependent dis- 
tributions appearing in Eq. (|8.7|) in the form 



etc. (8.11) 



K{x,PT)-f?(x,PT) 

= lft{x)-ff{x)]F?{x,p T ) 

+ f?(x) [F?(x,pt)-F»(x,pt)], 



If one integrated over px without restriction, as would 
correspond to integration of Eq. (|8.7[) over all Pt,tt, the 
difference of profile functions on the last line would in- 
tegrate to zero because of the normalization conditions 
Eq. (|8.10p . and the measured cross section difference 
would be proportional to the integrated valence quark 
density alone. However, without integration over all 
Pt the cross section difference contains an admixture of 
the antiquark distribution because generally F^{x,pt) 
F™(x,pt)- Consequently, one cannot extract the px~ 
integrated valence quark density from measurements of 
Eq. (|8.7j) over an incomplete range of Pt^ without ad- 
ditional assumptions about the intrinsic transverse mo- 
mentum dependence. This conclusion is general and not 
limited to the simple parton model used to illustrate the 
point. It also affects the use of the pion charge asymme- 
try (either as an absolute cross section difference or as a 
ratio of cross sections) in the flavor separation of polar- 
ized parton densities with proton and nuclear targets. In 
practice, there might still be considerable experimental 
advantages in using observables such as Eq. (|8.7[) for the 
study of pT-integrated parton densities; however, this re- 
quires detailed study based on the measured or modeled 
transverse momentum distributions. 

Generally, our findings underscore the importance of 
accurate measurements of the basic transverse momen- 
tum distributions of hadrons (tt, K) in unpolarized semi- 
inclusive DIS, up to Pr,h ~ 1 GeV, and differentially in 
x,z and Q 2 . Simple comparisons between the Pt dis- 
tribution of different particles (such as K + and K~) 
can serve as model-independent tests of the predicted 
pattern of intrinsic transverse momentum distributions. 
Measurements with a deuteron target would be partic- 
ularly useful, as they directly access the flavor-singlet 
(isoscalar) combination of the parton distributions and 
would allow one to search for signals of a difference 
between the transverse momentum distributions fi +d 
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and fi +d without the complications ol flavor separation. 
Such measurements provide essential information about 
the mechanism of particle production in semi-inclusive 
DIS and should be done before one studies more subtle 
observables such as spin asymmetries. 



C. Correlation measurements 

Much more information can be obtained from semi- 
inclusive experiments that measure particle production 
in the central and target fragmentation regions in corre- 
lation with an identified hadron in the current fragmenta- 
tion region (see Fig. ll9|) . Such measurements can answer 
the question of what "balances" the Pr,h of hadrons ob- 
served in the current fragmentation region, which is the 
key for unraveling the production mechanism in semi- 
inclusive DIS. As emphasized earlier, the Ptm of a hadron 
observed in the current fragmentation region can come 
from the parton intrinsic transverse momentum in the 
target, final-state interactions and QCD radiation in the 
hard process, or the soft fragmentation process. Single- 
inclusive measurements alone cannot discriminate be- 
tween these different sources. Correlation measurements 
offer additional observables that can test at least the rela- 
tive importance of the various mechanisms. A fully quan- 
titative theory of such measurements is a complex prob- 
lem and beyond the scope of the present article. Here we 
only wish to outline in what kinematic region such exper- 
iments could be performed such that they have a simple 
interpretation in terms of nonperturbative nucleoli struc- 
ture. 

In DIS experiments with measurements of multiparti- 
cle final states it is convenient to describe the produced 
particles in terms of their rapidity 
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1, Eh 
= — In 

2 E h 



P 



\\.h 



Pit 
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p2 
r T,h 



m" h y " is the hadron energy 



where Eh 

and P\\,h the longitudinal momentum, defined as the com 
ponent in the direction of the virtual photon momentum. 
The rapidity changes by a constant under Lorentz boosts 
along the longitudinal direction, so that rapidity differ- 
ences are frame-independent. The Lorentz-invariant ra- 
pidity interval over which particles with a given trans- 
verse momentum arc distributed is 



Y 



Vh,r 



Vh.i 



HW 2 /(P T . h +mt)}, (8.13) 



where W is the 7*7V center-of-mass (CM) energy and it 
is assumed that W 2 3> P T h + m\. Studies of hadronic 
final states in DIS at W ~ few GeV show that the target 
and current fragmentation regions occupy at least one 
unit of rapidity. To cleanly separate the two regions, 
rapidity intervals of Y ~ 4 are needed. Our picture 
based on QCD vacuum structure suggests that one look 
for correlations between pions with P% „ ~ 0.5 GeV 2 in 
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FIG. 19. (Color online) Measurements of particle production 
in the central and target fragmentation regions in correlation 
with an identified hadron in the current fragmentation region. 
y denotes the rapidity. 



the current fragmentation region and hadrons with sim- 
ilar transverse momenta in the target fragmentation re- 
gion. From Eq. (|8.13l) we see that this requires squared 
CM energies around W 2 « 30 GeV 2 . At these energies 
sea quarks with x = 0.05 (0.1) could be probed at scales 
Q 2 « xW 2 = 1.5 (3.0) GeV 2 , where perturbative QCD 
radiation still plays a minor role in inclusive particle pro- 
duction. The kinematic region described here thus rep- 
resents a "window" where one can expect to see nonper- 
turbative correlations between sea quarks reflected in the 
hadronic final state. 

By choosing a leading hadron in the current frag- 
mentation region with z > 0.5 (z is the fraction of 
the virtual photon's laboratory energy carried by the 
produced hadron) one can minimize also the effect of 
transverse momentum broadening through fragmenta- 
tion and approximately infer the intrinsic px of the struck 
(anti)quark in the target. By measuring particle pro- 
duction in the central and target fragmentation regions 
one can then establish how this transverse momentum is 
balanced in the final state. If the observed Pt,h of the 
forward hadron indeed originated from the nonperturba- 
tive intrinsic transverse momentum of quarks in the nu- 
cleon, the picture of pair correlations described in Sec. [VI] 
implies that the Pr,h should be balanced by the other 
(anti)quark in the pair, which materializes in the target 
fragmentation region, leading to long-range correlations 
in rapidity between the current and target fragmentation 
regions [79(. If, however, the main source of the observed 
Pt h were the fragmentation process, one expects short- 
range correlations within the current fragmentation re- 
gion. 

If Q 2 is increased from the values of ~ few GeV 2 dis- 
cussed above, QCD radiation should play an increasing 
role. In inclusive particle production this is described 
by DGLAP evolution and well understood theoretically 
and experimentally. In correlation measurements of the 
kind described here, these emissions should diminish the 
correlations expected from soft interactions by providing 
additional possibilities for the balancing of the observed 
Ptm in the current fragmentation region. 

Correlations between the current and target frag- 
mentation region could be studied by measuring either 
mesons or baryons, or a collection of hadrons, in the tar- 
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get fragmentation region. The role of baryons in balanc- 
ing the Pt,h of the current jet is likely more important 
in processes where a valence quark in the target is re- 
moved, potentially providing a handle to separate the 
intrinsic transverse momenta of valence and sea quarks. 

Correlation measurements of the kind described here 
would be feasible in the kinematic region covered by the 
CERN COMPASS experiment (squared lepton-nucleon 
CM energy s ~ 300 GeV 2 ) [80[ or with a medium-energy 
Electron-Ion Collider [8l|], if suitable detection capabili- 
ties for forward particles are provided. The merit of such 
correlation measurements in a smaller rapidity interval 
with the Jefferson Lab 12 GcV Upgrade (s tn 20GcV 2 ), 
especially with the CLAS12 detector, deserves further 
study |82| . 



D. Multiparton processes 

More direct tests of the idea of parton short-range cor- 
relations may become possible with the concept of multi- 
parton distributions, whose proper formulation in QCD is 
a subject of present work [70Ll7l|. Such distributions are 
required in the description of pp collisions with multiple 
hard processes, where correlations manifest themselves 
as an enhancement of the multiple process rate compared 
to the uncorrelated expectation. The data on three-jet 
plus phot on p roduction in pp collisions from the Teva- 
tron CDF (83| and DO [IH experiments show a significant 
enhancement compared to the uncorrelated expectation 
[HI , suggesting the presence of substantial correlations in 
the partonic wave function. Whether these correlations 
are due to perturbative QCD radiation or the nonpertur- 
bative vacuum structure discussed here is a challenging 
question which requires further study. 



E. Exclusive meson production 

The nonperturbative short-range correlations in the 
nucleon's partonic wave function discussed here (see 
Fig. I18p have potential implications beyond semi-in- 
clusive processes, for example for hard exclusive me- 
son production 7*iV — > M + N. In ir + production at 
Q 2 ~ few GeV 2 and x > 0.1, a substantial longitudinal 
cross section should result from the "knockout" of a cor- 
related qq pair in the nucleon with appropriate quantum 
numbers, which then forms a tt + by way of soft nonper- 
turbative interactions. This mechanism is theoretically 
sublcading compared to the hard scattering mechanism 
in the limit Q oo (in which the equivalent correla- 
tions are generated by perturbative gluon exchange) but 
may be practically dominant at all realistic values of Q 2 . 
It may also play a role in exclusive p° production in the 
region W ~ 2-4 GeV [86| , which shows an energy depen- 
dence consistent with the f-channel exchange of a qq pair 
with scalar quantum numbers [87| . The theoretical for- 
mulation of such knockout processes in exclusive produc- 



tion is the subject of on-going work. Generally, exclusive 
processes offer many possibilities for probing correlations 
in the partonic wave function due to QCD vacuum struc- 
ture, including their x-distribution and quantum num- 
bers. 
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Appendix A: Parametrization of soliton profile 

The profile function P(r) of the classical chiral field, 
Eq. (|2.6[) is determined by minimizing the classical en- 
ergy of the soliton, Eq. (|2.8p . Numerical minimization 
was performed with various types of cutoff of the Dirac 
sea contribution and result in very similar values of the 
classical energy; see Ref. [3l[ for a review. For a Pauli- 
Villars (PV) cutoff, described in detail in Sec. El the self- 
consistent profile was determined in Ref. [56[ . The result 
obtained with a dynamical quark mass M = 0.35 GeV, 
a PV cutoff with Mpy/M 2 = 2.52, cf. Eq. (j53^|) . and 
the chiral limit M v = 0, is shown in Fig. The figure 
shows the mesh points in the variable r at which the pro- 
file was determined by numerical minimization; at larger 
distances the profile effectively exhibits an asymptotic 
behavior as 

P{r) A/r 2 . (Al) 

With these parameters the classical energy of the soli- 
ton, i.e., the nucleon mass in leading order of the l/N c 
expansion, Eq. (|2.9p . was obtained as [56| 

M N w 3.26 M. (A2) 

The nucleon's isovector axial coupling, obtained by cal- 
culating the matrix clement of the axial current operator 
as a sum over quark s ingl e-particle levels, is gA = 1-03 
with these parameters [88| ■ In the chiral limit the isovec- 
tor axial coupling is proportional to the coefficient of the 
leading r~ 2 asymptotic behavior of the soliton profile at 
large distances, Eq. (|A1|) . 

9 A = y^ 2 A (A3) 
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FIG. 20. Soliton profile P(r) as function of the radius r. 
Points: Self-consistent profile obtained by minimizing the 
classical energy (PV regularization, M = 0.35 GeV, = 
0) [Hf|. Solid line: Parametrization P 8e if, Eqs. (|A4 J> — (| A7 p . 
Dashed line: Variational profile P var , Eq. f| A8f) . for M w = 0. 



The simple parametrization Eqs. (|A4[) - (jA7|) provides an 
excellent fit to the numerical self-consistent profile over 
all distances (see Fig. |20|). 

Also shown in the figure is the variational profile of 
Ref. [U, 

r 2 

P var (r) =-2arctan-|, (A8) 

r 2 = I.0A/- 2 , (A9) 

which was used extensively in earlier calculations of par- 
tonic structure [lU, [33j . This profile corresponds to 
the limit B — > oo of our more general parametrization, 
Eqs. (|A4[) and (|A5[) . and has an unphysical quadratic 
behavior near the center, which, however, is unimpor- 
tant for the quantities studied here (cf. the discussion in 
Sec. IV Cp . With the replacement 

2 2 

^ ^ (l + M„r)exp(-A/ ff r) (AI0) 

Eq. (|A8|) also allows one to study the effect of a finite 
pion mass. 

Appendix B: Bound— state level wave function 



and can alternatively be determined in this way. The 
equivalence of the two wa ys o f calculating gA for the self- 
consistent profile of Ref. [56| was verified in Ref. [§H . 

For the purposes of our study we need an analytic 
parametrization of the soliton profile which has the cor- 
rect limiting behavior at small and large distances and 
reproduces the numerical self-consistent profile with rea- 
sonable accuracy. It can be constructed in the form 



P S eif(r) 



-2 arctan T(r), 



(A4) 



where the function T(r) satisfies the following two condi- 
tions: (a) T(r) oc r _1 at r — > 0, to ensure linear behavior 
of the profile near the center; (b) r 2 T{r) = A + 0(r~ 4 ) 
for r — > oo, to guarantee the correct r~ 2 asymptotic be- 
havior of the profile and the absence of sublcading r~ A 
terms, which was noted in Ref. (30j in the context of a 
long-distance expansion of the equations of motion. A 
simple choice which fulfills these conditions is 



T(r) 



-| tanh(Br) 



(A5) 



The parameter r 2 , = A/ 2 determines the large-distance 
behavior of the profile; we fix it by requiring that the 
parametrization reproduce the "exact" numerical value 
of gA = 1-03 [88| via Eq. (|A3[) . which gives 



(A6) 



0.87 M- 



The parameter B is then fixed by fitting the behavior of 
the profile at small distances, which results in 



B 



1.6 M. 



(A7) 



In this appendix we give the explicit form of the 
bound-state level wave function and its Fourier trans- 
form, as used in the calculation of the level contribution 
to the transverse momentum distributions of quarks and 
antiquarks in Sec. IIVI Most of the relevant expressions 
were given in Appendix B of Ref. (32j |. and we include 
them here for reference only (see also Footnote [55| re- 
garding the convention for gamma matrices). 

In the standard representation of the gamma matrices, 
the bound-state solution of the Dirac equation in the 
nucleon rest frame, where the classical chiral field is given 
by Eq. (|2.6p . can be written in the form 



1 



V47T I - 



(xcr) 



h(r) 



|0>, 



(Bl) 



where r = \x\ and |0) is the spin-isospin wave function 
where spin and isospin arc coupled to zero total, 

fa + n)\o) = o, (0|0) = 1. (B2) 

The functions h and j are solutions of the radial equation 



d 2 

M cos P(r) - M sin P(r) 

or r 




The level wave function Eq. (|Bl[) is normalized to 



d 3 x $ 1 t ev (£c)$ lcv (a;) = 1, 



(B3) 



(B4) 
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corresponding to 



drr 2 [h 2 {r) + j 2 (r)] = 1. 



(B5) 



The wave function in momentum representation, de- 
fined according to Eq. (|3.29[) , 

*iev(p) = Jd 3 xe- l P x $i ev (x), (B6) 
can be written in the form 

« / %) \ 

$ lev (p) = ^ (p<r) .. . |0), (B7) 

where p = \p\. The radial functions in the momentum 
representation here are given by 



h{p) = / drr 2 h(r)R p0 (r), 



j(p) = J drr j(r)R pl (r), 
'() 



(B8) 



where R p i are the free radial wave functions of the con- 
tinuous spectrum, defined as 

R pl (r) = J^Ph(pr) = ^/IdVshw (B9) 
V 7r V 7r V P dr J r 

where = 0, 1, . . .) denote the spherical Bessel func- 
tions, and normalized according to 



drr 2 R pl (r)R p , l (r)=5(p-p'). (BIO) 



The normalization condition for the momentum repre- 
sentation of the bound-state wave function, Eq. (|B7[) . is 



d 3 p 



*lev(p) t $lov(p) = 1, 



corresponding to 



dp[h 2 (p)+j 2 (p)] =1. 



(Bll) 



(B12) 



When computing matrix elements between bound- 
state level wave functions, one can take advantage of the 
properties of the spin-isospin singlet state. Specifically, 



.r'|0) = ...(-a i )|0), 

|r\.. = (0|Hx i; )-.., 
l|cr l |0) = (0|r*|0) = 0. 



(B13) 
(B14) 
(B15) 



Using these identities one can convert matrix elements of 
products of a and r matrices into equivalent expressions 
involving only one kind of matrix, which can then be 
evaluated using standard techniques. 

Appendix C: Fourier transform of soliton field 

When evaluating the sea quark distribution using the 
gradient expansion we need explicit expressions of the 
3-dimensional Fourier transform of the static classical 
chiral field in the nucleon rest frame. In coordinate space 
the field is given by Eq. (|2.6|) and can be expanded as 



U c \{x) = cosP(r) + - — ^sinP(r), 



(CI) 



where r = \x\. The Fourier transform Eq. (|5.7|) is of the 
form 

U cl (k) = Jd 3 x er lkx [U cl (x) - 1] 



47T 



8(k) + { -^ P (k) 



(C2) 



where k = |fe|. The scalar functions s(k) and p(k) are 
determined as 

/>oo 

s(k)= / drr 2 j (kr) [cos P(r) - 1] , (C3) 



p(k) 



drr 2 ji(kr) sinP(r), 



(C4) 



where jo an d ji are the spherical Bessel functions, cf. 
Eq. dSSJ. The traces in Eq. ([021 and Eq. (pTTTj) are 
then obtained as 

tr fl [f7 cl (fe)C/ cl (fe)t] = 327r 2 [.s 2 (fc) + p 2 (k)}, (C5) 

fc 3 



tr fl [r 3 i7 cl (fc)C/ cl (fe)t] = 64tt 2 s(k) p{k) 



k 



(C6) 



In proving the large-iV c inequalities for the sea quark 
distributions, cf. Eq. (|5.77[) . it is important to note that 
the sum and difference of the forms in Eqs. (|C5|) and (|C6[) 
is non-negative. Namely, 

tr fl [(l±T 3 )[/ cl (fc)[/ cl (fc)t] 



= 32tt 2 
> 32?r 2 



! (fc) + p 2 (k) ± 2s{k)p(k) — 



\ s (k)\ 2 + \p(k)\ 2 -2\ S (k)\\p(k)\^-l 



Ifc 3 ! 



>32^ 2 (| S (fc)|-b(fc)|) 2 , 
where the last inequality follows from |fc 3 | < k. 



(C7) 
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